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Abstract 

We consider pure three-dimensional quantum gravity with a negative cosmological 
constant. The sum of known contributions to the partition function from classical geome- 
tries can be computed exactly, including quantum corrections. However, the result is not 
physically sensible, and if the model does exist, there are some additional contributions. 
One possibility is that the theory may have long strings and a continuous spectrum. An- 
other possibility is that complex geometries need to be included, possibly leading to a 
holomorphically factorized partition function. We analyze the subleading corrections to 
the Bekenstein-Hawking entropy and show that these can be correctly reproduced in such 
a holomorphically factorized theory. We also consider the Hawking-Page phase transition 
between a thermal gas and a black hole and show that it is a phase transition of Lee- Yang 
type, associated with a condensation of zeros in the complex temperature plane. Finally, 
we analyze pure three-dimensional supergravity, with similar results. 
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1. Introduction 

This paper is devoted to describing some explicit computations relevant to three- 
dimensional pure quantum gravity with negative cosmological constant. The classical 
action can be written 

Three-dimensional gravity with a cosmological constant was first studied in [[[]. Some 
early milestones in the study of this theory, both of them special to the case of negative 
cosmological constant, have been the construction of an asymptotic Virasoro symmetry 
and the recognition that the theory admits black holes The asymptotic Virasoro 
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symmetry is now understood as part of the structure of a dual two-dimensional conformal 
field theory || and this structure together with modular invariance can be used to deter- 
mine the entropy of a black hole of asymptotically large mass |J. For a recent review of 



None of the results mentioned in the last paragraph are special to the case of pure 
three-dimensional gravity, without additional fields. Our intent here, however, is to study 
this minimal theory, with the goal of computing its exact energy levels in a spacetime that 
is asymptotic at spatial infinity to the classical "vacuum" of three-dimensional Anti de 
Sitter space, AdS3. As always in General Relativity, a well-posed problem is obtained by 
specifying what the world should look like at spatial infinity - in this case we ask that it 
should coincide with AdS3 at infinity - and then analyzing all possible "interiors." Our 
problem in this paper is to compute the precise quantum energy levels that arise. We 
perform a computation based on known concepts, and in a sense that we will explain, this 
computation is not successful. The reasons for this failure are not clear and we consider 
several hypotheses. 

It is convenient to summarize the spectrum of energy levels in the form of a trace 
Tr exp(—/3H), where H is the Hamiltonian and f3 is a positive real number (or more 
generally a complex number with positive real part). As usual in General Relativity, the 
Hamiltonian is defined via the ADM procedure in terms of the leading behavior at spatial 
infinity of the correction to the pure AdS3 metric. It was indeed by carefully examining 
this procedure that the asymptotic Virasoro symmetry (which generalizes the obvious 
conservation laws such as conservation of energy) was discovered . 

There is also a conserved angular momentum J which generates a rotation at infinity 
of the asymptotically AdS3 spacetime and commutes with H. Consequently, one can 
introduce an additional parameter 9 and try to compute a more general partition function: 



This partition function is naturally computed via a Euclidean path integral. Accord- 
ing to the standard recipe, the integral is taken over Euclidean three-geometries that are 
conformal at infinity to a two-torus £ with modular parameter r = 6/2-k + if3. We write 
\dz\ 2 for a flat metric on this torus, where z is a complex coordinate subject to the identifi- 
cations z — *■ 2 + 1, z — > z + t, and introduce another coordinate u > 0, such that conformal 
infinity will be the region u — > 0. Then the metric should behave for u — > like 



the extensive work on this subject, with references, see 0. 



Z(p, 6) = Tt exp(-(3H - id J). 



(1.2) 



ds 2 



dz\ 2 +du 2 
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plus subleading terms. 

The Euclidean path integral also has a natural interpretation in the dual two- 
dimensional CFT [||: it is the genus one partition function of this dual theory, on the 
two-dimensional surface E. In fact, by definition, the Hamiltonian and momentum H and 
P of the dual CFT coincide with H and J, the Hamiltonian and angular momentum of the 
original theory in AdS3. This dual interpretation is informative, but is not really needed 
to motivate the computation that we perform. 

This Euclidean path integral is a formal recipe, for various reasons. One problem is 
that, in general, the Euclidean quantum gravity path integral is not convergent because 
the action is not bounded below ||. The only known way to deal with this problem is 
to expand around a classical solution; in doing so, one can obtain a meaningful result at 
least in perturbation theory. There is no clearly established claim in the literature that 
topologies that do not admit classical solutions do not contribute to the Euclidean path 
integral. But there is also no known method to evaluate their contributions. 

It turns out that in the present case, one can describe the classical solutions completely. 
This is explained in Sec. 2.1. The classical solutions are precisely the ones considered in 
U and in |I(J (in studying the elliptic genus of certain AdS3 models derived from string 
theory). Moreover, perturbation theory around a classical solution terminates with the 
one-loop term, and that term can be easily evaluated by adapting the arguments of @. 
This is the content of Sec. 2.2. 

We are therefore in a position to write down the complete sum of known contributions 
to the path integral. This is done in Sec. 3. The sum turns out to require some regular- 
ization, and we use what we believe is a natural regularization, analogous to zeta function 



regularization ||11|| . We ultimately obtain an explicit, though complicated, formula for the 
sum of known contributions to Z(f3, 9). However - and this is our main result - the sum 
is not physically sensible: it cannot be written as Tr exp(—f3H — i9J) for any commuting 



operators H, J in a Hilbert space. (According to jl2;1 , this possibility was also conjectured 
by S. Minwalla.) 

We do not know the correct interpretation of this result. In Sec. 4, we discuss some 
possibilities. One possibility is that the minimal theory that we postulate here actually 
does not exist. There are many quantum theories that do exist that look semiclassically 
like three-dimensional gravity coupled to additional fields; indeed, there are a plethora of 



known string-derived models, such as the ones studied in [10]. However, it is not clearS 
that there exists nonperturbatively a minimal theory that one would want to call pure 
gravity. One possible interpretation of our result is to indicate that this minimal theory 
does not exist. 

The other possibility, broadly speaking, is that there are additional contributions 
to the path integral beyond the known ones that we evaluate in Sees. 2 and 3. To 
be concrete, we consider in Sec. 4 two scenarios. One is that the minimal theory of 
three-dimensional gravity, in addition to the known BTZ black holes and Brown- Henneaux 
boundary excitations, also describes cosmic strings. The motivation for this proposal is 
that known models of three-dimensional quantum gravity, such as the string-based models 



considered in |I(J , do always have cosmic strings. Perhaps this is also true for the "minimal" 
theory, if it exists. 

The second scenario is that in addition to the real saddle points that are classified and 
evaluated in Sees. 2 and 3, an exact description of the theory should also include complex 
saddle points. We describe a specific scenario in which the inclusion of complex saddle 
points leads to the holomorphically factorized partition function that was proposed (based 
on highly speculative reasoning) in [[13 j . This involves a doubled sum over saddle points 



similar to what is assumed in [|12| , |r4}| . The resulting partition function is consistent with 
an interpretation as Tr exp(— f3H — i6J) with some Hilbert space operators H, J. 

Though we are not able to put this proposal in a convincing form, we do uncover 
one interesting fact, which concerns the semiclassical behavior of the partition function 
(assumed to be holomorphically factorized and extremal) as the gravitational coupling G 
goes to zero with fixed AdS radius I. In this semiclassical limit, the partition function 
Z = Tr exp(—(3H — i9J) is always dominated, as long as (3 and 9 are real, by a real saddle 
point corresponding to a real geometry. The partition function can be dominated by a 
classical saddle point, but only if one asks a more exotic question with complex values of 
(3 and 6. 

In Sees. 5 and 6 we discuss questions of black hole physics in three dimensions. In 
section 5 we discuss black hole entropy: a theory of pure quantum gravity, assuming it 



1 The problem may be analogous to trying to define a minimal string theory in four dimensions. 
There are many theories that macroscopically are four-dimensional string theories, such as theories 
obtained by compactification to four dimensions from the critical dimension, or gauge theories 
with flux tubes or vortex lines, but it is plausible that there is nothing that should sensibly be 
called a minimal four-dimensional quantum string theory. 
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exists, should give a proper microscopic accounting of the Bekenstein-Hawking entropy of 
the BTZ black hole. In the semi-classical limit this entropy is just the horizon area. The 
computations of section 2, however, allow us to determine the perturbative corrections to 
this semi-classical result, which are described in section 5.1. In section 5.2 we show that 
an infinite series of corrections can be reproduced in a holomorphically factorized theory. 

In investigating these holomorphically factorized theories, we noticed an additional 
interesting phenomenon that is the subject of Sec. 6. Three-dimensional quantum gravity, 
as a function of 0, exhibits in the semi-classical limit the Hawking-Page phase transition 
|15|| between a thermal gas (in this case, a gas of Brown-Henneaux boundary excitations) 
and a black hole. At first sight, it is puzzling how such a phase transition can be compatible 
with holomorphic factorization. We show, however, that this question has a natural answer, 
in terms of a condensation on the phase boundary of Lee- Yang zeroes of the partition 
function. 

Finally, in Sec. 7, we extend our analysis to pure Af = 1 supergravity, obtaining 
results similar to what we find in the bosonic case. 

2. Known Contributions To The Path Integral 

Our task in this section is two-fold. The first step is to classify Euclidean solutions of 
Einstein's equations, with negative cosmological constant, with the asymptotic behavior 
described in the introduction. We write M for the three-dimensional spacetime, and £ for 
its conformal boundary; as in the introduction, E is a Riemann surface of genus 1. We 
assume that M is smooth, that the metric on M is complete, and that £ is the only "end" 
of M. According to the standard logic of Euclidean quantum gravity, known contributions 
to the path integral have these properties; it may be appropriate physically to relax them, 
but we do not know how. 

We proceed in two steps. In Sec. 2.1, we describe the possible choices of M. The 
result is actually well-known in the theory of hyperbolic three-manifolds. In Sec. 2.2, 
we evaluate the contribution to the partition function of a particular M. The sum over 
different choices of M is postponed to Sec. 3. In Sec. 2.3, we describe the general form of 
the partition function in a theory with finite entropy. 
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2.1. Classification Of Solutions 

The automorphism group of AdS3 is >SO(3, 1), which is the same as SL(2, C)/Z2. We 
may write the metric on a dense open subset of AdS3 asi 

ds 2 = \ dz \ 2 + du \ w>0 , zeC. (2.1) 
u 

If we combine the (z, u) coordinates into a single quaternion y = z + ju, the action of an 
element (* G SX(2,C) can be written succinctly as 

y^{ay + b){cy + d)- X . (2.2) 

In this expression the element -i) e SL(2,C) acts trivially, so ( |2.2|) actually de- 
scribes the action of SL(2, C)/7*2 on AdS3. In general, any classical solution M of three- 
dimensional gravity with negative cosmological constant takes the form AdS 3 /r, where T 
is a discrete subgroup of SO(3, 1) and AdS 3 is the part of AdS3 on which T acts discretely. 

The conformal boundary of M can be constructed as follows. First of all, the conformal 
boundary of AdS3 is a two-sphere, which one can think of as CP 1 , acted on by SL(2, C) in 
the usual way. This CP 1 may be regarded as the complex z-plane in ([O]) at u — > plus 
a point at infinity. From (|2.2|), one can see that SL(2, C) acts on this CP 1 in the familiar 



fashion 

z -> — . (2.3) 

cz + a 

To construct the conformal boundary of M = AdS 3 /r, one first throws away a certain 
subset of CP 1 in a neighborhood of which the discrete group T acts badly. This set is closed, 
so its complement is an open subset U C CP 1 . 

The discrete group V acts freely on U, and the conformal boundary of M is the 
quotient 

E = U/T. (2.4) 

Since the action of SX(2,C) preserves the holomorphic structure of CP 1 , U/T carries a 
natural holomorphic structure and the isomorphism ( |2.4|) is valid holomorphically, not 
just topologically. 

Given this, let us investigate the condition for E to be of genus 1. E is topologically a 
two-torus, so its fundamental group is 7Ti(E) = Z © Z. It follows from eqn. fl2.4p that the 



2 We set £ = 1 unless otherwise indicated. 
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fundamental group of U is a subgroup of the fundamental group of E. Possible subgroups 
tti(U) C Z x Z are of three types: 

(i) tt\{U) may be a subgroup of 7Ti(E) of finite index, isomorphic to Z©Z. (A special 
case is iri(U) = 7Ti(E).) 

(ii) tvi(U) may consist of multiples of a single non-zero vector x G 7Ti(E), in which 
case 7Ti(J7) = Z. 

(mj 7Ti(C7) may be trivial. 

In case (i), U is a finite cover of E, and therefore is itself a Riemann surface of genus 
1. However, a Riemann surface of genus 1 is not isomorphic to an open subset of CP 1 . So 
case (i) cannot arise. 

Cusp Geometry 

In case (Hi), U is the universal cover of E, and so is isomorphic to C (or P 2 ). The 
holomorphic structure of C is unique up to isomorphism. C is isomorphic holomorphically 
to an open subset of CP 1 in essentially only one way: it is the complement of one point, 
say the point at z = oo. The subgroup of SL(2, C) that leaves fixed the point at infinity 
consists of the triangular matrices 

The point zGC corresponds to ^ ^ ^ G CP 1 , and a triangular matrix acts by z — > X 2 z+Xw. 

Since U is simply-connected, V must be isomorphic to Z © Z (in order to get the right 
fundamental group for E). Any discrete group of triangular matrices that is isomorphic to 
Z © Z is generated by two strictly triangular matrices 

{I ;)■ (J !)■ (-> 

where a and 6 are complex numbers that are linearly independent over R. (If a and b 
are linearly dependent over R, then the group generated by these matrices is not discrete, 
since a suitable linear combination ma + nb, for m, n G Z, can be arbitrarily small.) Up 
to conjugacy by a diagonal matrix, the only invariant of such a group is the ratio b/a. 
Therefore, we can reduce to the case a = 1, b = r, and moreover by taking b — > —6 (which 
does not affect the group generated by the two matrices) we can assume that Imr > 0. 
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We have therefore arrived precisely at the group of symmetries 



z — > z + m + nr, m, n G Z (2-7) 

of the complex z-plane. The quotient is a genus 1 surface E with an arbitrary r-parameter. 
The three- manifold M = AdS^/T can be described very concretely and is a standard 
example. The (u, z) coordinates in (|2.1| ) cover precisely the subspace AdSg on which T 
acts "nicely." So M = AdS^/T is given by the metric (|2.1|), now subject to the identification 

But M does not obey the conditions described at the beginning of this section. M is 
smooth and has a complete Einstein metric. However, in addition to the "end" at u = 0, 
which is the one required by our boundary condition, M also has a second "end" at u = oo. 
Our problem was to classify Einstein manifolds with only one end; M does not qualify. 

The second end of M is one at which the metric of E collapses to zero (rather than 
blowing up, as it does for u — > 0). An end of this kind is known as a "cusp." In the theory 
of hyperbolic three-manifolds, one often considers hyperbolic metrics with such cusps. One 
of the main results about them is that, under suitable conditions, such a metric can be 
slightly perturbed so as to eliminate the cusp. In the present context, this gives the sort 
of metric that we will describe momentarily. 

Physically, we cannot be certain that omitting the case of the "cusp" is the correct 
thing to do. However, restricting to cuspless metrics appears to be the right thing to 
do in known examples of AdS3 theories. And pragmatically, we believe that from what 
is known of three-dimensional gravity, it would be difficult to give a sensible procedure 
for evaluating the contribution of the spacetime with the cusp. The reason for the last 
statement is that non-trivial one-cycles (loops in E) become sub-Planckian in length near 
a cusp, so a semiclassical treatment is not valid. 

Semiclassical Geometries 

The last case to consider - case (ii) - is that the fundamental group of U is Z. This 
means that topologically U is K x S 1 (S 1 is a circle). The holomorphic structure of U is 
then uniquely determined to be that of the z-plane minus a point, which we may as well 
take to be the point at z = 0. The subgroup of triangular matrices that preserve the point 
z = is simply the group of diagonal matrices. 
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Therefore T is a discrete subgroup of the group of diagonal matrices. It cannot be a 
finite group (or the first Betti number of £ = U/T would be 1, while the desired value is 
2). This being so, there are essentially two cases to consider. 

First, r may be isomorphic to Z, generated by a matrix of the form 

W= (l A J ESL(2X)- (2.8) 



q- 

By exchanging the two eigenvalues we can assume that \q\ < 1. (If |g| = 1, then either q is 
a root of unity and the group generated by If is a finite group, or q is not a root of unity 
and the subgroup of SL(2, C) generated by W is not discrete.) Alternatively, F may be 
isomorphic to Z x Z n , generated by W together with 

f exp(27rz/n) \ , . 

\ exp(-27rz/n) J ' 1 ' 

with some integer n. 

Let us first consider the case that T = Z. Then E = U/T is obtained from the complex 
z-plane by throwing away the point z = and dividing by the group generated by W. It 
is convenient to write z = exp(27riu>), so that w is defined modulo 

+ 1 (2.10) 

and W acts by 

«; ^ «; + — — . (2.11) 
2m 



The quotient of the ty-plane by (|2.10|) and (|2.11|) is a Riemann surface of genus 1, as 
required. The complex modulus of this surface is r = -y^, i.e. it is given by q = e 2 ' KlT . 
More generally, however, the modulus of this Riemann surface is defined only up to r — > 
(or + b)/(cr + d) with integers a,b,c,d obeying ad — be = 1. Therefore, we will get an 
equivalent Riemann surface if 

q = exp(2yrz(ar + 6)/(cr + d)) (2.12) 

for such a, 6, c, d. 

One might conclude, therefore, that we get a three-manifold obeying the required 
conditions for every choice of a, 6, c, d. This is not quite the case, for two reasons. First, an 
overall sign change of a, 6, c, d does not affect q or the associated three-manifold. Second, 
once c and d are given, a and b are uniquely determined by ad — be = 1 up to shifts of 



9 



the form (a, b) — > (a, 6) + £(c, d), £ 6 Z. Under this transformation, g as defined in eqn. 
( |2.12|) is invariant. So the possible three- manifolds really only depend on the choice of the 
pair c, d of relatively prime integers, up to sign. For each such pair, we find integers a, b 
such that ad — be = 1 , and identify q via eqn. (|2.12|) . This gives a manifold that we will 



call M c d- This family of manifolds were first discussed in the context of three-dimensional 
gravity in ||. 

The Geometry of M c ^ 




Figure 1: a) An infinite cylinder representing AdS3. The boundary of the 
cylinder represents conformal infinity; time translations act by vertical shifts, 
b) A slice of height determined by (3. The manifold M^i is built by gluing 
together the top and bottom, after a rotation that identifies the boundary 
points marked by solid dots. 

The simplest such manifold is M^i, which we will now describe in more detail. If we 
identify the real one-parameter subgroup diag(e 6 ,e _b ) of SX(2,C) as the group of time 
translations, the AdS3 metric ( |2.1|) can be put in the form 

ds 2 = cosh 2 r dt 2 + dr 2 + sinh 2 r d<f> 2 , (2.13) 

with — oo < t < oo, < r < oo, and < (p < 2n. Described in this coordinate system 
is the subset AdS' 3 C AdS3 on which F acts nicely; its topology is D x R, where D is a 
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two-dimensional open disc parameterized by r and <p. Conformal infinity is at r = oo. (See 
Fig. 1.) The element diag(e b , e~ b ) acts by t — > t + b. The group of spatial rotations is the 
one-parameter group diag(e l61 , e~ l9 ), acting by </> — > </> + 6- 

The group element W therefore generates a combined time-translation and spatial 
rotation. To explicitly divide by W, we "cut" AdSg at times t = and t = 2nlmr. Then 
we glue together the top and bottom of the region < t < 2n Im r after making a spatial 
rotation by an angle 27rRer. This is sketched in Fig. lb. The resulting spacetime Mo i is 
topologically D x S 1 , so its fundamental group is indeed Z. 

The path integral in this spacetime has a simple semiclassical meaning, since it may 
be interpreted in terms of Hamiltonian time evolution. A state is prepared at time zero 
and propagates a distance f3 = 2nlmr forward in Euclidean time. In this process, the 
state vector is multiplied by the time evolution operator exp(— [3H), where H is the Hamil- 
tonian. Then, after a spatial rotation by an angle 9 = 2n Re r, which acts on the state 
by exp(— i#J), we glue the top and bottom of the figure, which results in taking the in- 
ner product of the final state with the initial state. The whole operation gives the trace 
Tr exp(—/3H — i9J) defined in the Hilbert space of perturbative fluctuations around AdS3. 
This Hamiltonian interpretation of the path integral of Mq^ will be the basis for evaluating 
it in Sec. 2.2. 

The other manifolds M c ^ are obtained from M^i by modular transformations, that 
is, by diffeomorphisms that act non-trivially on the homology of S. This fact will allow us 
to evaluate their contributions to the path integral. These other manifolds may be thought 
of as Euclidean black holes. 

For example, in the (t, <p) coordinates introduced above, Mq^ involves the identifica- 
tions 

+ it ~ + it + 2tt ~ (p + it + 2tvt. (2.14) 

From the metric (|2.13 ), note that the <p circle is contractible in the bulk, since the coefficient 
of d<p 2 becomes zero at the origin r = 0. The manifold M^o is described by the same 
coordinate system (|2.13|) , but with new identifications 

2tt 

(p + it ~ (p + it + 2n ~ (p + it . (2.15) 

r 

These identifications ( |2 . 1 5| ) may be written in the form of ( p. 14] ) by taking (j) + it — > 
- (<fi + it). After this transformation, the coordinate that is contractible in the bulk now 
involves a combination of d> and t. 
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If we take r to be imaginary, the scaling by 1/r has the effect of exchanging (and 
rescaling) <fi and t. is hence a Euclidean black hole, in which the "time" circle is 

contractible rather than the "space" circle. If we rotate to Lorentzian signature by t — > it, 
then the locus r = where the coefficient of dt 2 vanishes is the horizon of the black hole. 



In fact, Mi ; o is just the Euclidean version of the BTZ black hole [T^]. The more general 



manifolds M c ^ are often referred to as the SL(2, Z) family of black holes. 
Orbifolds 

Finally, let us consider the extension in which F has a second generator given by 
fl2.9|) . We may assume that n = 2m is even, since the element diag(— 1, — 1) (the non- 
trivial element of the center of £X(2,C)) acts trivially on CP 1 . We also may as well 
assume that m > 1, since the case m = 1 leads to nothing new. 

What we get when V has an additional generator with m > 1 is simply a three- 
dimensional space of the form M Cj d/Z m . The conformal boundary is still a Riemann 
surface of genus 1, and by changing q we can adjust its modular parameter as we wish. 
However, for m > 1, the group element ( |2.9| ) acts on AdS' 3 with fixed points, meaning 
that M Cjt //Z m has orbifold singularities. The fixed points are of codimension 2 and the 
singularities look locally like 1R 2 /Z m , where Z m acts as a rotation by an angle 2n /m. This 
produces a deficit angle 

6 = 2tt(1- 1/m). (2.16) 

The picture is sketched in Fig. 2. 

At the classical level, the physical meaning in three-dimensional gravity of a codimen- 
sion two singularity characterized by a deficit angle is usually that it represents the orbit 



of a massive particle. The mass of the particle is related to the deficit angle [fETf. The 
particular values of the deficit angle in eqn. ( |2.16|) are special because they are related to 
orbifolds. 

Of course, there are many consistent theories of three-dimensional quantum gravity 
plus matter, obtained from various widely studied string theory and M-theory construc- 
tions. They each have massive particles of various sorts, the precise masses and spins being 
model-dependent. It would be hopeless to try to completely solve such a general theory. 

The question of interest in the present paper is whether there exists a theory that is 
in some sense minimal and solvable. Our hypothesis is that the minimal theory should be 
described by smooth geometries without orbifold or deficit angle singularities. Pragmati- 
cally, it appears difficult to avoid the problem found in Sec. 3 simply by including a small 
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Figure 2: The black line represents an orbifold singularity in the interior 
of a spacetime that at infinity looks like AdS 3 . 

set of such singularities. In any event, were we to allow singularities, we would not know 
which ones to allow. 

In the dual CFT, massive particles correspond to primary operators of positive dimen- 
sion. Any two-dimensional CFT - and they are abundant, of course - can be interpreted as 
a dual AdS3 theory. If this dual theory has a macroscopic, semiclassical interpretation, it 
can be interpreted in terms of AdS3 with particles and black holes. From the holographic 
point of view, it is difficult to understand which CFTs have a semiclassical interpretation, 
but it is clear that, whether or not there is what one might call a minimal theory, there 
are many consistent theories with massive particles. 

2.2. Evaluation Of The Partition Function 

We write Z c ^{r) for the contribution to the partition function of the manifold M Cj d- 
Because the manifolds M Cy a are all diffeomorphic to each other, the functions Z Cj d(r) can 
all be expressed in terms of any one of them, say Zo 5 i(r), by a modular transformation. 
The formula is simply 

Z c , d (r) = Z 0A ((ar + b)/(cr + d)), (2.17) 
13 



where a and b are any integers such that ad — be = 1. The partition function, or rather 
the sum of known contributions to it, is 



Z(r) = Z cAr) = E Z ^ aT + 6 )/( cr + ( 2 - 18 ) 

c,d c,d 

The summation here is over all integers c and d which are relatively prime and have c > 0. 

This formula shows that the key point is to evaluate Zq^(t). We recall that this 
contribution is simply Tr exp(—(3H — i9J), computed in the Hilbert space that describes 
small fluctuations about AdS3 (as opposed to black holes). If we know the eigenvalues of 
the commuting operators H and J in the Hilbert space of small fluctuations, then we can 
compute the trace. 

In the most naive semiclassical approximation, Zq \(t) is just exp(— J), where I is 
the classical action. In computing this action, one can not just evaluate the action fll.lQ 
for the solution ( |2.13| ); such a computation would give an infinite answer, coming from 
the boundary at r — > oo. The full action includes the Gibbons-Hawking boundary term, 
which has the opposite sign of the Einstein-Hilbert term ( |1 . 1| ) . This extra term removes 
the divergence, and one arrives at a finite (negative) answer for the action of M^i | fL6f : 

I = -47vklmr (2.19) 

where k = £/16G. Therefore, in this approximation, we have 

Z 0> i(r) ~ \qq\~ k . (2.20) 

The result ( |2.2U| ) has a simple interpretation, in terms of facts explained in . Three- 



dimensional pure gravity with the Einstein-Hilbert action ( |1 . 1| ) is dual to a conformal field 
theory with central charge cj_, = cr = 31/ 2G = 24k. The formula cl = cr = 31/ 2G is 
actually the semiclassical approximation of @ . Depending on how the theory is regularized, 
there may be quantum corrections to this formula, but they preserve cl = cr because 
the Einstein-Hilbert theory is parity symmetric. (By adding to the action a gravitational 
Chern-Simons term, we could generalize to cl ^ cr. We omit this here.) For our purposes, 
we simply parametrize the theory in terms of k = c^/24 = cr/24. Since cl and cr are 
physical observables (defined in terms of the two-point function of the stress tensor in 
the boundary CFT, or an equivalent bulk computation), the theory when parametrized 
in this way does not depend on any choice of formalism. (By contrast, if we describe 
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the results in terms of the microscopic variable £/G, then it is also necessary to describe 



the regularization.) It is conjectured in |13[] that pure three-dimensional quantum gravity 
exists, if at all, only for integer k. However, for our present purposes, we do not need to 
know if this is correct. 

Let Lq and Lq be the Hamiltonians for left- and right-moving modes of the CFT. 
They are related to what we have called H and J by 

H — L + L 

(2-21) 

J — Lq — Lq. 

The CFT ground state has Lq = — cl/24, Lq = — cr/24, or in the present context 
Lq = Lq = —k. Equivalently, this state has H = —2k, J = 0. Its contribution to 
Tr exp (— 27r(Im r)H + 27rz(Rer) J) is exp(47r/clmr) = \qq\~ h , as in eqn. ( 2.20 ). 



Thus, although seemingly only a naive approximation, eqn. (|2.20|) actually gives the 
exact result for the contribution of the ground state to the partition function (if the theory 
is parametrized in terms of the central charges). As a result, it is also exact for Imr — > oo. 

To get a complete answer, we also need to know the energies of excited states. In more 
than 2 + 1 dimensions, the trace would receive contributions from a gas of gravitons as well 
as from other particles, if present. Their energies would in general receive complicated 
perturbative (and perhaps nonperturbative) corrections. Hence, although the function 
Zq j i(t) has a natural analog in any dimension, above 2 + 1 dimensions one would not 
expect to be able to compute it precisely. 

In 2+1 dimensions, there are no gravitational waves. Naively, there are no perturbative 
excitations at all above the AdS3 vacuum, and hence one might at first expect the formula 
(|2.20|) to be exact. However, this is mistaken, because of the insight of Brown and Henneaux 
0. There must at least be states that correspond to Virasoro descendants of the identity, 
or in other words states obtained by repeatedly acting on the CFT vacuum |0) with the 
stress tensor. If L n and L n are the left- and right-moving modes of the Virasoro algebra, 
then a general such state is 

oo oo 
n=2 m=2 

with non- negative integers u n ,v m . (For a state of finite energy, almost all of these integers 
must vanish. The products begin with n, m = 2 since L_i and L_i annihilate the CFT 
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ground state.) A state of this form is an eigenstate of L and L with L = — k + ^2 n=2 nu n , 
Lq = —k + J2m=2 mv m- The contribution of these states to the partition function is then 



>o,i(r) = l^r fc nr _ 2| j_ gn|2 - (2-23) 
It is convenient to introduce the Dedekind r\ function, defined by 

oo 

V(r)=q 1 ^l[(l-q n )- (2-24) 



n=l 



Eqn. ( J2.23| ) can then be rewritten 

^o,i(r) = ^|^|-( fc - 1 / 24 )|l- ? | 2 . (2.25) 

This formula will be useful in Sec. 3 because (Imr) 1 / 2 |r7(r)| 2 is modular-invariant. 

We will show shortly how the Virasoro descendants of eqn. ( |2.22|) arise in an anal- 
ysis along the lines of Brown and Henneaux. The argument will also show that the for- 
mula ( 2.23Q is exact to all order of perturbation theory. As for whether there are non- 
perturbative corrections to eqn. (|2.23|) , it is not clear whether this question is well-defined, 
since it may be impossible to separate the question of unknown nonperturbative correc- 
tions to Zq^i from the more general question of unknown nonperturbative corrections to 
the exact partition function Z. At any rate, it will be easier to discuss what nonpertur- 
bative corrections would mean after explaining why eqn. (|2.23| ) agrees with perturbation 
theory. 

Before presenting this argument, let us discuss what the answer means from the general 
point of view of quantum mechanical perturbation theory. Usually, one defines an effective 
action I e ff such that the partition function is Z = exp(— /c/ e ff )• -^eff is equal to the classical 
action / plus contributions generated from r-loop diagrams: 

oo 

Jeff =I + J2 k ~ ri r- (2-26) 

r=l 

Here I r is a function generated from Feynman diagrams with r loops. In a general theory, 
the expansion (|2.26[ ) is only asymptotic in 1/k; there may be further contributions that are 
non-perturbatively small for k — > oo. However, in the present context, the formula ( |2.23| ) 
implies that (for perturbation theory around the manifold Mo,i) there are no corrections 
beyond one-loop order; I e g is simply the classical action plus a one-loop correction. In other 
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words, the formula implies that if Zq^ is directly computed in perturbation theory, then 
the perturbation series terminates with the one- loop term. (Such a direct computation is 
also briefly discussed below.) 

In a sense, this should not come as a surprise. The gauge theory interpretation of 



three-dimensional gravity [|T8|JT9[1 (or even more naively the absence of any local excitations) 
suggests that in some sense it is an integrable system. It is often true for quantum integrable 
systems that properly chosen quantities are one-loop exact. 

Derivation Of The Formula 

In most of physics, the quantum Hilbert space of a theory is constructed by quantizing 
an appropriate classical phase space. This framework may ultimately be inadequate for 
quantum gravity, but it is certainly adequate for constructing the perturbative Hilbert 
space that is needed to compute Z 0j i in perturbation theory. 

The phase space of any physical theory is simply the space of its classical solutions. 
In 2 + 1-dimensional gravity with AdS boundary conditions at infinity, naively speaking 
the phase space relevant to perturbation theory consists of only a single point, since any 
solution that is close to the AdS3 solution (close enough to exclude black holes) actually 
is diffeomorphic to AdS3. 

However, we must be more careful here. In General Relativity, one should divide only 
by diffeomorphisms that approach the identity fast enough at infinity. After doing so, one 
constructs the classical phase space Ai; it parametrizes classical solutions that obey the 
boundary conditions modulo diffeomorphisms that vanish fast enough at infinity. One is 
then left with an action of a group G that consists of those diffeomorphisms that preserve 
the boundary conditions, modulo those that vanish fast enough at infinity that they are 
required to act trivially on physical states. 

In General Relativity in 3+1 or more dimensions, in a spacetime that is asymptotically 
Lorentzian, the group G is the Poincare group. This is why quantum gravity, in an 
asymptotically Lorentzian spacetime, is Poincare invariant. In 2 + 1-dimensional gravity 
with negative cosmo logical constant, the obvious analog of this answer would be SL(2, R) x 
SX(2,R) (or a group locally isomorphic to that one), the group of symmetries of AdS 3. 
However, Brown and Henneaux showed that the actual answer turns out to be the 
infinite-dimensional group G = diff S 1 x diff S 1 (which contains <SX(2,R) x SX(2,M) as a 
subgroup). By computing Poisson brackets, they also showed that the quantum symmetry 
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group will be not G but a central extension with (in a semiclassical approximation) cl = 
cr = 3£/2G. We denote this central extension as diff S 1 x diff S 1 . 

We will carry the analysis of [|J just slightly farther to describe the phase space Ai 
and the resulting energy levels. First of all, A4 is a homogeneous space for G, since if 
we divide by all diffeomorphisms, then the classical solutions that obey the boundary 
conditions (and are close enough to AdSs) are equivalent to AdS3. So M. is G/H, where 
H is some subgroup of G. In fact, H is the subgroup of G that leaves fixed a given point 
on M.. Differently put, the symmetry group of a point on M. is isomorphic to H. In 
the present context, we simply pick the point on M. corresponding to AdS3 and observe 
that its symmetry group is SX(2,R) x SL(2, M). So this is H, and the phase space is 
M = (diff S 1 x diff 5 1 )/(5L(2, E) x SL(2,R)) = (diff S 1 /SL(2, R)) 2 . 

In general, the quantization of a homogeneous space G/H should be expected to give 
a Hilbert space that is an irreducible representation R of G (or more generally of a central 
extension of G). Moreover, as G/H contains an i7-invariant point, the representation R 
will contain a vector that is an eigenvector for the action of H (and hence is .ff-invariant 
if if is a simple non-abelian group). 

In the present problem, a representation of the Virasoro group diff S 1 in which Lq is 
bounded below and in which there is an SL{2, R)-invariant vector is uniquely determined 
once the central charge c = 24k is given. It is the "vacuum" representation, containing a 
vector |0) with (L n + kd n fi)\Q) = 0, n > — 1. This representation is spanned by states of 
the form 

oo 

UL%\n), (2.27) 



n=2 



with energy 



e = -k + ^2na n . (2.28) 

n=2 

Quantizing (diff S 1 /SL(2, R)) 2 gives a tensor product of two such representations for the 
two factors of G — (diff S 1 ) 2 , and in this way we arrive at the spectrum claimed in eqn. 
( 2.22|) , and hence at the formula ( |2.23| ) for the partition function. 

Quantization of homogeneous spaces of diff S 1 is described in [ 20[| . The fact that 
quantization of the quotient diff 5' 1 /5'L(2,R) gives the vacuum representation of diff S 1 is 
explained from several points of view in eqns. (170), (172), and (174) of that paper. 

The states that we have described can be regarded as boundary excitations, supported 
near the boundary of AdS3, since M. collapses to a point if we are allowed to make general 
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changes of coordinate near the boundary. We call these boundary excitations the Brown- 
Henneaux (BH) states, since they are so closely related to the analysis in 0. 



Corrections ? 

Now let us ask to what extent there may be corrections to the spectrum just described. 

First we consider the question of quantum corrections to the formula (|2.28|) for this 
family of states. We can dispose of this question immediately: this particular family 
of states transforms in an irreducible representation of the symmetry group, and this 
representation has no possible deformations once the central charge is fixed. 

So there are no quantum corrections to the energies of the states just described. Could 
there be additional states contributing to Zq^I 

In general, in three-dimensional gravity, there may be additional states of various 
kinds, corresponding for instance to massive or even massless particles. What do we expect 
in "pure" three-dimensional gravity, if it exists? It is important to distinguish two length 
scales: the AdS3 length scale £, and the Planck length G. In the semiclassical regime of 
k — > oo, we have G << I. 

In eqn. ( |2.28[ ), energies are measured in units of l/£. A reasonable minimum condition 
for anything that one might call a theory of pure gravity is that the energy of any state 
that is not part of the BH spectrum goes to infinity in the semiclassical limit k —>■ oo (when 
measured relative to the energy of the ground state). States that are not part of the BH 
spectrum and whose excitation energies remain fixed for k — > oo would be interpreted in 
terms of non-gravitational fields propagating in AdS3. 

As for what sort of states might have energies that go to infinity for k —>■ oo, we 
certainly expect BTZ black holes, with energy bounded below by k. As for what else there 
may be, the sky is the limit in terms of conceivable speculations. If so inclined, one can 
postulate solitons with excitation energies proportional to k, states similar to D-branes 
with energies proportional to /c 1 / 2 , etc. Orbifold singularities with fixed deficit angle, as 
described in relation to eqn. ( |2.16| ), would also have excitation energy of order k. 

Our point of view is that the question of whether there are nonperturbative corrections 
to Zq^i from states whose excitation energy diverges as k —>■ oo is hard to separate from 
the more general question of unknown nonperturbative corrections to the exact partition 
function Z. Since the perturbative evaluation of Zq^ gives a convergent and physically 
sensible result (which is even one- loop exact), we may as well regard hypothetical contribu- 
tions in which Mo 1 is enriched with a soliton, a D-brane, an orbifold singularity, or some 
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other unknown type of excitation as representing different, presently unknown sectors of 
the path integral. In this sense, the formula for Zo y i is exact. 



Comparison To Perturbation Theory 

Since we claim that the formula for Zq^ is one-loop exact, the question arises of why 
not to simply calculate it by evaluating the relevant one- loop determinants. 



We know of two efforts to do so. In [16|, a formula much simpler than our result 



( p.23|) is claimed. The analysis relies on the relation of the relevant product of determi- 



nants to Ray-Singer and Reidemeister torsion [JTJ. For an irreducible flat connection on a 
compact manifold without boundary, the torsion is simply a number, but for a manifold 
with boundary (such as M 0i i effectively is), the torsion must be understood as measure on 
a certain moduli space associated with the boundary. A rather subtle treatment of this is 
needed, we suspect, to compute the one-loop correction using its interpretation via torsion. 

On the other hand, in [FZT], the one-loop correction was studied by expressing the 
determinants in terms of the appropriate heat kernels, which were evaluated using a method 
of images (analogous to the derivation of the Selberg trace formula) . We believe that this 
method is conceptually completely correct. The claimed result is qualitatively similar to 
our formula (p.23|) (for example, it has an expansion in integer powers of q and q) , and we 
hope that it will prove possible to resolve any discrepancies between the result in [^T| and 
the one claimed here. 

As for why it is much simpler to compute the one-loop correction via the Hamiltonian 
route that we have followed, this should not really come as a surprise. In general, path 
integrals on a product S 1 xY are often most easily evaluated by constructing an appropriate 
Hilbert space in quantization on Y and then taking a trace.i That is especially likely to be 
true in the present situation, in which the relevant excitations are rather subtle boundary 
excitations, already known but difficult to rediscover. 



3 An excellent example of this statement is given by Chern- Simons gauge theory with a compact 
gauge group for the case that Y is a compact Riemann surface without boundary. The partition 
function on S 1 x Y is an integer, the dimension of the physical Hilbert space 7i associated with 
Y . It is not too hard to describe TC, compute its dimension, and thereby learn the value of the 
path integral on S 1 x Y. But direct evaluation of this path integral by Lorentz-covariant methods 
is difficult, even in perturbation theory. 
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2.3. The General Form Of The Partition Function 

We conclude with an analysis of the following question: in general, what sort of 
function can be written as Tr exp(—(3H), where if is a hermitian operator on a Hilbert 
space HI We assume that H is constrained so that Tr exp(— (3H) is convergent whenever 
Re (3 > 0. (According to standard assumptions about pure three-dimensional gravity, this 
condition is satisfied, but we discuss in Sec. 4.1 one way that it might fail.) 

For Tr exp(—(3H) to be convergent, the number he of states of energy less than E 
must be finite for each E. Indeed, the trace is bounded below by % exp(— (3E). Finiteness 
of he for all E implies in particular that H must have a discrete spectrum. 

Let E* be any value of the energy and let Ex, . . . , E n be the energy eigenvalues that 
are no greater than E*. Their contribution to the partition function is J27=i ex P(~ P^i)i 
and the full partition function is therefore 

n 

Tr exp(-pH) = ^exp(-/3£,) + 0(exp(-PE*)). (2.29) 

i=l 

Here C(exp(— fiE*)) is a function that is bounded by a multiple of exp(— (3E*). 

We have to be careful here with one point: any finite sum of states of energy greater 
than E* makes a contribution of the form 0(exp(— /3E*)), but in general an infinite set 
of states with energy > E* may make a contribution that is not bounded in this way. 
However, our hypothesis that Tr exp(— (3H) converges whenever Re/3 > ensures that the 
number he of states of energy less than E grows with E more slowly than any exponential. 
This is enough to justify the error estimate in eqn. ( |2.29| ). 

Our problem is slightly more general. We have a pair of commuting hermitian op- 
erators H and J, and we are computing the trace Tr exp(—(3H — iOJ). However, the 
eigenvalues of J are integers, so 

Tr exp(-/3# - i6J) = ^e~ %je Tr^- exp{-/3H), (2.30) 

where Tr, is a trace in the subspace Tij in which J acts with eigenvalue j. The functions 
Trj exp(— j3H) should be constrained exactly as in eqn. ( |2.29| ). 

In Sec. 3, we will find instead that the evaluation of known contributions to 
Tr exp(— (3H) take the form just described up to a certain energy E*, beyond which 
this form breaks down. The result looks like 

n oo 

^expi-PEi) + eM-PE*)(f + J2 ( 2 - 31 ) 
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where / is a negative integer, rather than being positive as it should be, and the corrections 



decay as power laws, rather than exponentials. (The series in ( |2.31|) may be only asymp- 
totic.) Both the fact that / is negative and the fact that there are power-law corrections 
mean that this function cannot be interpreted as Tr exp(— (3H). 

Finally, let us compare the assertion that Tr exp(— f3H) can only converge if H has a 
discrete spectrum with standard results about the thermodynamics of physically sensible 
systems with continuous spectrum, such as a free gas of particles on the real line. What is 
usually computed in such a case is the free energy per unit volume; the total free energy is 
infinite simply because the volume is infinite. To make the total free energy finite, one can 
place the gas in a finite volume; then the spectrum is discrete and Tr exp(— j3H) converges. 
A gas of particles in Anti-de Sitter space will have a discrete spectrum and finite partition 
function Tr exp(— f3H) as well. This is because the motion of any finite energy particle in 
Anti-de Sitter space is restricted to a finite volume region in the interior of AdS. 



3. Computing the Sum over Geometries 

As explained in Sees. 1 and 2, the known contributions to the partition function 
of pure gravity in a spacetime asymptotic to AdS3 come from smooth geometries M Cj d, 
where c and d are a pair of relatively prime integers (with a pair c, d identified with 
— c, —d). Their contribution to the partition function, including the contribution from the 
Brown-Henneaux excitations, is 



Z(t) = X>,i(7t), 



(3.1) 



where 



and 



7T 



ar + b 
ct + d' 



Q 



n=2 



7 



a b 
c d 



e SL(2,Z) 



q' 



m\ 



\v(r) 



(3.2) 



(3.3) 



The summation in ( |3.1| ) is over all relatively prime c and d with c > 0. Since Zo t i(r 
is invariant under r 
and b in 



r + 1, the summand in (|3.1| ) is independent of the choice of a 
This sum over c and d in (|3.1| ) should be thought of as a sum over the 
coset PSL(2, Z)/Z, where Z is the subgroup of PSL(2, Z) = SL(2, Z)/{±1} that acts by 
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r — > r + n, n G Z. Given any function of r, such as Zo i i(r), that is invariant under 



r — > r + 1, one may form a sum such as (|3.1[) , known as a Poincare series. 



The function v / Imr|r7(r)| 2 is modular-invariant. We can therefore write Z(t) as a 
much simpler-looking Poincare series, 

Z{r) = * £ (VWI^I-^ 1 / 24 !! - ,| 2 ) , (3.4) 

where (. . .) | is the transform of an expression (. . .) by 7. Writing out explicitly |1 — q\ 2 = 
1 — q — q + qq, we see that we really need a sum of four Poincare series, each of the form 

E(t; n, m) = (V^q~ n q~ m ) , (3.5) 

c,d 

with n — m equal to or ±1. Precisely this sum, or rather its s-dependent generalization 
introduced below, has been studied in Sec. 3.4 of as was pointed out by P. Sarnak.S 
K we set k = n + m, \i = m — n, and use the fact that Im (77*) = Imr/|cr + <i| 2 , then the 
basic Poincare series can be written 

E(t; k, h) = Vlmr |cr + <i| _1 exp {2itk Im 7T + 27rz/x Re 7r} . (3.6) 

c,d 

When k = and \x = 0, this sum is a non-holomorphic Eiseinstein series of weight 1/2. 
Sometimes we omit r and write just E(k,[j,). In terms of this function, the partition 
function is 

Z(t) = * (£(2/c - 1/12, 0) + E(2k + 2 - 1/12, 0) 

VlmTl^T)! 2 (3.7) 

- E(2k + 1 - 1/12, 1) - E(2k + 1 - 1/12, -1)). 

5.1. TTie Regularized Sum 



The sums ( |3.6|) that define the Poincare series we need are divergent. (Such diver- 
gences have been encountered before in similar sums related to three-dimensional gravity 

ig^3H0 Usin § 



ar + b _ a _ 1 ^ 
ct + d c c(cr + d) ' 



Actually, the generalization to arbitrary integer values of n — m is considered in |22]. We 



would encounter the same generalization if we modify the original Einstein-Hilbert action (1.1) 
to include the gravitational Chern-Simons term. 
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one may show that if r = x + iy then 

y 



Im (77-) 



[cx + d) 2 + c 2 y 2 
a cx + d \ • ) 

Re ( 7 r) = - - c (( cx + d y + C 2 y 2y 

For \i = 0, the exponential factor in the definition of E is 

exp {2ttk Im 77-} = exp \ 2ttk -| — ^ I , (3.10) 

^ (cx + d) 1 + c l y l J 



and this goes to 1 for c, d — > 00. So the first two terms in ( |3.7| ) diverge linearly as 
S c d l cr + a ^ l ar § e c an d ^- The other terms also diverge, though more slowly. 

We claim that this divergence has a natural regularization. On the upper half plane, 
which we call IK, there is a natural SL ( 2, R) -invariant Laplacian: 

A = -!/ 2 f|USV (3.11) 



dx 2 dy* 

A short calculation shows that the function y 1 ! 2 is an eigenfunction of A: A(y 1 ^ 2 ) = 
(1/4)?/ 1 / 2 . Since A is £X(2, IR)-invariant, the same is true of (Im.'jT) 1 ^ 2 for any 7 G SL(2, Z) 
(or even SX(2,R)): 

A a/Iiu (77-) = iv/lm(7r). (3.12) 

Using this, one may verify that although the Poincare series for E(r; k, fj,) is divergent, 
the corresponding series for (A — 1/4)E(t; k, h) actually converges. (This series is obtained 
by acting termwise with (A — 1/4) on the Poincare series for E.) To see this, we just 
observe, using ( pj.lUp , that contributions in which derivatives appearing in A act on the 
exponential factor in ( |3.6| ) get an extra convergence factor of order l/(max(c, d)) 2 and lead 
to a convergent sum over c and d. On the other hand, contributions in which none of the 
derivatives act on the exponential actually vanish, because of ( |3.12j ). 



So (A — 1/4)E requires no regularization. Similarly, if we set F = Vim r | rj \ 2 , then 
(A — 1/4) (FZ) requires no regularization, where Z is the partition function. Since no 
regularization is required, we assume that the naive sum gives (A — 1/4) (FZ) correctly. 

The Laplacian A, acting on the Hilbert space of square-integrable £X(2, Z)-invariant 
functions on !K, has a continuous spectrum starting at 1/4. It also has a discrete spectrum, 
but there are no discrete modes with eigenvalue 1/4. So the operator A — 1/4 is invertible 
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acting on square-integrable functions, and roughly speaking, we now want to argue that 
once (A — 1/4) (FZ) is known, FZ is uniquely determined. 

There is a subtlety in making this argument, because actually Z and FZ are not 
themselves square-integrable. Rather, we expect that Z grows exponentially for y — > oo, 
because the AdS3 vacuum has negative energy — 2k. There are also Brown-Henneaux 
excitations of this vacuum, again with negative energy, and they contribute additional 
exponentially growing terms for y —>■ oo. If we assume that all exponentially growing 
contributions to Z come from the known Brown-Henneaux states, then this, together with 
a knowledge of (A — 1/4) (FZ), is enough to determine Z uniquely. 

But for our main application, we do not really need this assumption (which appears to 
be not quite correct according to the analysis of Sec. 3.3). The explicit calculation of Sees. 
3.2 and 3.3 will give us a function Z, compatible with the Poincare series for (A — 1/4) (FZ), 
such that the fastest growing exponentials agree with the Brown-Henneaux spectrum. The 
leading departure from the contributions of that spectrum are determined at the end of 
Sec. 3.3 and take the form 

1 / (7r 3 -67r)(ll + 24fc) , 

~ 9C(3l V 

5(53tt 6 - 882tt 2 ) + 528(tt 6 - 90tt 2 )A; + 576(tt 6 - 90tt 2 )/c 2 _ 2 _ 3 , 

+ 2430C(5) V (2/ ' 

(3.13) 

This function cannot be written as a sum of exponentials with positive integer coefficients, 
since the leading coefficient is negative, and the corrections do not have the right form. 
Any other candidate for Z would be obtained by adding a correction Z such that (A — 
1/4) (FZ) = 0. But that equation is not obeyed by any function Z that would solve our 
problem; the equation (A — 1/4) (FZ) = is not obeyed by —Z' or by any function that 
differs from it for y — > oo by a sum of exponentials with positive integer coefficients. So 
the problem that we will find is not affected by adding to the function Z that we compute 
an additional contribution Z that obeys (A — 1/4) (FZ) = 0. 

More Convenient Alternative 

Although this argument is satisfying conceptually, it does not give a convenient way 
to determine Z in practice. A much more convenient method is to adapt ^-function 



regularization [II] to this problem. 
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In the present context, the analog of (^-function regularization is to replace the Poincare 
series ( |3.5|) by a more general one depending on a parameter s: 

E(t; s, n,m) = J2 ((Imr) V n g _m ) | 7 ■ (3.14) 

c,d 

Since (Im7r) s = y s /\ct + d\ 2s , the series E(t; s, n, m) converges for Res > 1. Our original 
problem concerns the case s = 1/2. As we will see shortly, E(t; s, n, m), defined initially 
for Res > 1, can be analytically continued to s = 1/2 without any problem. This analytic 
continuation gives a natural way to define the original function E{r; n, m) and therefore 
the partition function Z. This turns out to be a very practical and useful method to study 
E(r; n, m). 

The only problem with this approach is that the physical meaning of the parameter 
s in three-dimensional gravity is unclear; hence, it is not clear a priori that the analytic 
continuation will give the right answer. The argument involving (A — 1/4) (FZ) is clearer 
conceptually, because (A — 1/4) (FZ) is a physical observable and we simply use the fact 
that the formal path integral expression for it converges. 

However, it is not difficult to show that the two methods give the same result. For 
this, we observe that the function H(r; s, m, n) = (A + s(s — l))E(r; s, m, n) can also be 
represented by a Poincare series, convergent when Res > 0. For Res > 1, where both 
Poincare series converge, the two functions E and H, both defined by their Poincare series, 
obey H = (A + s(s — 1))E. This automatically remains true after analytic continuation 
to s = 1/2. So the function E(t; 1/2, m, n) defined by analytic continuation from Re s > 1 
has the property that (A — 1/A)E(t; 1/2, m, n) is given by the obvious Poincare series. 

In terms of k = m + n, (jl = m — n, the regularized Poincare series is 

E(s, k, a) = I l9 exp {2nn Im 7T + 27ri// Re 77-} . (3.15) 

\CT + d\ Zs 

c,d 

When k = and \i = 0, this sum is the non-holomorphic Eisenstein series of weight s. 



3.2. Poisson Resummation 



Sec. 3.4 of |22| contains precisely what we need to analyze the sum ( |3.15| ), make the 
analytic continuation, and determine if the partition function Z is physically sensible. 
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First, define d = d' + nc, where n is an integer and d! runs from and c — f . We may 
separate out the sum over n in ( |3.15| ) to get 



E(s, = y' e ^y+i^) + E/( c ' d » ( 3 - 16 ) 

c>0 d'el/cZ nel 



where 



, y , 2,71 Ky fa cx + d 

f(c, d,n)= exp + 2tt^ - 



c(r + n) + d'\ 2s \ |c(r + n) + d'\ 2 \c c\c{r + n) + d'\ 2 

(3-17) 

The first term in eqn. ( 3.f 6 ) comes from c = 0, d = 1. 

The Poisson summation formula allows us to turn the sum over n into a sum over a 
Fourier conjugate variable n 

^/M>) = ^;(c jC f,n) (3.18) 
where f(c,d',n) is the Fourier transform 



/oo 
dn e 2 ™/(c,d',n) 
-oo 



- -p (« - -)) /_: * u^)*- 

(3.19) 

We have written the integral in terms of a shifted integration variable t = n + x + — . Upon 
Taylor expanding the exponential that appears in the integral and introducing T = t/y, 
we get 

;m» = £ e ^rt^-) ^ jy e - (« - 

m=0 '"' 

(3.20) 

In ( |3.2U| ), c and d! do not appear in the integrals but only in the elementary prefactors, 
so we can study the sums over c and d' explicitly. Note that, as fi is an integer, ( |3.20| ) 
depends on a only modulo c. The value of a mod c is determined by d', given that ad' = 1 
mod c. For a given d', such an a exists if and only if d! lies in the set (Z/cZ)* of residue 
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classes mod c that are invertible multiplicatively. So, dropping the prime from d' , we may 
write the sum over that variable as 

S(-h,fJL;c)= J2 ewl^i( ~ nd+ ^ d - )}, (3.21) 

where d~ x G (Z/cZ)* is the multiplicative inverse of d. This sum is known as a Kloosterman 
sum. 

Rearranging the sums in ( |3.16| ), we now have 



E(s, k, fi) = y s e 2n{Ky+lflx) + e'^^Eftis, k, y) (3.22) 

n 

where 

oo / oo \ 

E*(s,k,h) = W(*>«>/*) y x - m - s \ Y^c- 2 ^ + ^S{-h^-c) . (3.23) 

m=0 \c=l / 

Here we have defined the integral 



ml J-oo 



(3.24) 



Note that ( |3.23|) is independent of x, so that (|3.22|) has the form of a Fourier expansion 



in x with Fourier coefficients Ea(s, k, fi) given by (|3.23| ). These Fourier coefficients are 
typically complicated functions of y, since the integral ( |3.24j ) depends on y. 



These Fourier coefficients are precisely what we want. In view of ( |3.71 ), the Fourier 
expansion of the function E with respect to x will give a similar expansion for the partition 
function Z. The Fourier coefficients of Z are the functions Tr-, exp(—(3H), the partition 
function restricted to states with angular momentum J = j. These are the functions that 
we want to understand. 

We are almost ready to analyze what happens when we continue the formulas to 
s = 1/2. For m > 0, the integral in (|3.24|) is convergent for Res > 0, and likewise the sum 



J2c~ 2im+s) S(-h^;c) (3.25) 



c=l 



converges for Re s > 0. A problem does occur for m = 0, since then neither the integral 
nor the sum is convergent at s = 1/2. We return to this shortly. 
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3.3. Theh = Mode 



Let us first consider the Fourier mode which is constant in x, i.e. the h = term in 
( p.22|) . In this case the integral ( |3.24| ) is independent of y, and may be evaluated explicitly. 
For /i = 0, the result can be expressed in terms of V functions 

„ 2 m 7r m+1 / 2 r( s + m - 1/2) 

lm,o(a,«,0 = n m / 3.26 

m\T(s + m) 

while for jj, = ±1, we require also hypergeometric functions 

J m ,o(-, «, ±D - cos (-) ^Tsf (— +*> 7 + ' ; 2^ 

TO7 r\ (2tt)"T (f ) T (f + s) (1-m m 3 2 



+ m.sin — j w , r(w + g) 2F1 , -j + s;-; k 

(3.27) 

When m = 0, this formula simplifies to 

/ ,o(^^±l) = v^ r(g r ~^ /2) - (3-28) 



The sum over c may also be evaluated exactly (see Sec. 2.5 of [p22j). For /x = this 
evaluation involves the Kloosterman sum 5(0, 0; c); from the definition ( |3.21[ ), one can see 
that S(0, 0; c) is equal to the Euler totient function 0(c) (which is defined as the number of 
positive integers less than c that are relatively prime to c). The sum over c is a standard 
one 

± „-*»+• >S(0, 0; c) = £ c- 2 <»*+^(c) = C( ^ ( ( ;' ( 'J;; ) - ) 1) . (3.29) 

This formula can be obtained as follows. We start by noting a basic property of the totient 
function: for any n, ^2^\ n <p(d) = n. In order to evaluate the sum ^2 c c~ a 0(c), let us 
multiply this sum by ((a) = ^ n~° ' . This gives 



oo 

(nc) 

c=l n,c=l 
oo 



C(a)^c-V(c)= J> 



J^m 47 J] 0(c) (3.30) 

m=l c|m 



oo 



= X)"» 1 - <F = C(^-i). 



m=l 
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Setting a = 2(s + m) gives (|3.29| ). 

For [x = ±1, the Kloosterman sum becomes a special case of what is known as Ra- 
manujan's sum (see Sec. 2.5 of P2"|): 

S(n,0;c)= Yl e 2nid/c = fi(c). (3.31) 

de{i/ciy 

Here n(c) is the Mobius function, which is defined as follows: (i[c) = if c is not square- 
free, while if c is the product of k distinct prime numbers, then (i(c) = (— l) k . The sum 
over c is given by 

oo oo 1 

X>- 2 ^S(0,±l;c) = X>- 2 ^V(c) = c(2(g + m)) - (3.32) 

To prove this, one may use the basic property of the Mobius function: for any n, 
Sd|n M(^) = &n,\- To compute X] c c_cr A i ( c )' we mu hiply this sum by £(er) to get 

oo oo 

C(a)^c->(c)= ^(cn)-V(c) 

c=l n,c=l 
oo 

= J2 m ~ a Y,v( c ) (3-33) 

"1=1 c |m 

oc 

= m" cr 5 r71j i = 1. 

m=l 

Setting a = 2 (to + s) gives (|3.32|) . 

Putting these formulae together gives an exact expression for the x independent part 
of our Poincare series ( |3.16| ) : 

oo 

E (s, K,n)= 22 w m(s, n, n)y l ~ m ~ s . (3.34) 

m=0 

The constants w m (s, k, fi) in this expansion are independent of x and y, and are given 
explicitly by equations ( |3.26| ), ( |3.27| ), (|3.29| ) and ( |3.32| ) - we will write them out in more 
detail for the cases of interest below. The expansion ( |3.34j ) is one of our key results. It is 
an explicit series expansion of our Poincare series in powers of y~ x 

Now we can study the behavior at s = 1/2 of the delicate contribution with m = 0. 
For (i = 0, this term is finite because the factor of T(s — 1/2) appearing in ( |3.26| ) and 
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the factor of ((2s) appearing in (|3.29|) both have simple poles at s = 1/2; as s — > 1/2, 
r(s — l/2)/C(2s) — > 2. In fact, all of the complicated factors cancel to give 

E Q (l/2,n,0) = -y 1 / 2 + O(y- 1 / 2 ). (3-35) 

The coefficient of the leading term is negative because £(0) = —1/2. For \i = ±1, the 
m = term is finite because the factors of V(s — 1/2) in Q3.28|) and C(2s) in ( 3.32 ) both 
have simple poles at s = 1/2. We find 

#,(1/2, K, ±1) = 2y 1/2 + 0(y" 1/2 )- (3-36) 

Going back to (|3.22|) , and writing simply E(k, fi) for £'(1/2, «, //), we now have 

E(k, 0) = y 1/2 exp(27my) - y 1 / 2 + 0(y~ 1/2 ). (3.37) 

Actually, we have not yet discussed the Fourier modes with n^O; however, the integrals 
( |3.24| ) vanish exponentially for y — > oo and so do not affect the assertion in eqn. ( |3.37| ). 
Likewise, we have 

E(k, ±1) = y 1/2 exp(2yr(Ky ± ix)) + 2y 1/2 + 0(y~ 1/2 ). (3.38) 

If we evaluate eqn. ( |3.7| ) for the partition function Z keeping only the first terms in 
eqns. ( |3.37| ) and ( |3.38| ) - the exponentially growing terms - then the formula for Z simply 
reduces to Zq^, the contribution of the Brown- Henneaux states. What we have gained 
from all the work that we have done is that we can now calculate corrections to Zo,i- 
The leading corrections come from the corrections to the exponential terms in ( |3.37| ) and 
( |3.38| ). Adding them up and evaluating (|37?|), we get 

Z = Z ,i + (-6 + Oiy- 1 )) . (3.39) 

Since 1 

/M 2 ~ kl" 1/12 , z is governed by the Brown-Henneaux spectrum up to energy 
—1/12 (slightly below the classical black hole threshold, which is at zero energy). The 
number of states at that energy is not a positive integer, as one would hope, but rather 
-6. 

Moreover, the derivation shows that the corrections are given by a power series in 1/y, 
not a sum of exponentials. As was explained in Sec. 2.3, for Z to have an interpretation 
as Tr exp(—f3H), the corrections must be given by exponentials. 
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Thus we have arrived at the main conclusion of this paper: the sum of known con- 
tributions to the partition function of pure three-dimensional gravity is not physically 
sensible. 

It is illustrative to compute the next few terms in the expansion of the partition 
function (|3.39| ) in powers of y. This means taking m ^ in equations (|3.26| ), ( |3.27|) , ( |3.29| ) 
and ( |3.32| ). Unlike the m = case, these expressions do not have poles at s = 1/2, so 
the computation is straightforward. Let us start with the \x = case. We find that the 
coefficients w m (s, K, fx) are given by 

/ / s 2 m 7r m+1 /V(2m) 
w m 1/2, «, = - TT^ ■ 3.40 

ml (to + l/2)Q{2m + 1) 

So the next few terms in the Poincare series (|3.34|) are 



/ 2tt 3 \ _ 1/9 . / 4 



E (l/2, K , 0) = -y 1/2 + ™ K y -i/a + ^ y -3/ 2 + 0(y -5/ 2) (3 _ 41) 



Uc(3) K r" + Vl35C(5) 

To evaluate the [i = ±1 terms, note that at s = 1/2 the first two arguments of the relevant 
hypergeometric function appearing in ( |3.27 ) are integers. This means that the formula 
([3.27 ) simplifies considerably at s = 1/2 - it is just a polynomial in k. It is 

1 27T m+1 /2 

I m>0 -, «, ±1) = , . T m (K) (3.42) 

2 ml (to + 1/2) 

where T m (k) denotes a Chebyshev polynomial of the first kind. So the coefficients appear- 
ing in ( p.34| ) are 

w m (l/2, «, ±1/2) = 7^ r m (»c). (3.43) 

ml (m + l/2)(\2m + 1) 

This allows us to write down the next few terms in the series: 
£ (l/2, k, ±1) = 2y 1/2 + ( 2T 1/2 + ( ^(2^ 2 - 1)1 2T 3/2 + 0(y" 5/2 ). (3-44) 



.c(3) ; y vsas) 

Going back to (|3.22| ), and writing E(k, fi) for E(l/2, k, /a), we have 

e(k,o) = y ^eM^^y)~y 1/2 +{^f) y~ 1/2 +(ji^f 2 ) y- 3/2 +o(y~ 5/2 ) (3-45) 



and 



E(k, ±1) = yV*exp(2n(Ky ± zx))+2y 1 /2 + y -i/a + ^_^ (2 ^ - i)) y -3/ 2+0(y -6/ 2) 

(3.46) 
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Plugging this into (|3.7|) gives the expansion of the partition function 

7-7 + 1 ( fi I (^ 3 -M(ll + 24fc) 

Z - Zo > 1 + W^~ 9C(3) V 

5(537r 6 - 8827r 2 ) + 528(7r 6 - 907r 2 )/c + 576(7r 6 - 907r 2 )/c 2 _ 2 ^. 

H — — v + 0(y~ 

2430C(5) y {y 

(3.47) 

The additional contributions to the partition function in this expression have two notable 
features. First, and most importantly, they are not zero. Thus, as described above, the 
partition function truly cannot be represented as a sum of exponentials. Second, they 
differ qualitatively from the leading y° term: the additional coefficients appearing here are 
positive and irrational, rather than negative and integer. 

34. h^O Modes 



Now we will consider the h ^ terms. For \i = and h ^ 0, the integral ( pj.24j ) is a 
-KT-Bessel function 

I m>fl (s, k, 0) = 1 J y s+m - 1 / 2 K s+m _ 1/2 (27r|n|y). (3.48) 

mil (s + m) ' 

The Kloosterman sum is now the general case of Ramanujan's sum 

S(n,0;c)= e 2 ™ d / c = J2^)- (3-49) 

d^(Z/cZ)* S\n 

We will simply quote the answer for the sum over c (see Sec. 2.5 of [f22"lQ 

^ C(2(* + m))^ 

Taking s = 1/2 gives a Fourier coefficient of the regularized partition function. Con- 
sider first the m = term. For n = 0, this was the dangerous term in the analytic 



continuation, but for h 7^ 0, it simply vanishes, because the C(2s) in (|3.50|) has a pole at 
c = 1/2, causing the Kloosterman sum to vanish, and - unlike the ft = case - the integral 
( |3.48| ) is finite at s = 1/2. The other terms are non-zero, and give 

00 2 3 / 2 7 r rn + 1 lf?l m / \ 

E h (l/2,k,0) = if? .4-1/9W9 Z^n E 5 " 2m ^m(27r|n|y). (3.51) 

^ m!r(m + l/2)C(2m + 1) I ^ / 
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Each term in the sum over m vanishes for large y as e _27r ' n ' y . 

Unfortunately, when // ^ the coefficients in the expansion of the partition function 
are more difficult to compute. The problem is not the integral ( |3.24j ); these integrals are 
finite and can be evaluated analytically, although we will not write the answer here. The 



answer is a finite sum of Bessel functions of the form appearing in (|3.48|) , each of which 
is multiplied by a polynomial in y. As y — > oo, these Bessel functions cause the integrals 
vanish as e _27r l n l y , just as in the \i — case described above. 
However, when \x 7^ 0, the sum over c 

00 

£ c -a(™+'>S(fi, M , c ), (3.52) 

c=l 

though of considerable number-theoretic interest, cannot be expressed in terms of familiar 
number-theoretic functions such as the Riemann zeta function. Analytic properties of 
these sums have been extensively studied. We will simply quote the relevant results, ft 
has been shown that the sum ( |3.52| ) defines a meromorphic function on the complex s 
plane. This function is essentially the Selberg zeta function associated to the modular 
domain D = 9C/SL(2,Z). When we take s = 1/2, the function ( 3.52 ) remains regular for 
elementary reasons if m > 0. Indeed, one can see directly that the sum converges for these 



values, by noting from the definition of the Kloosterman sum ( 3.21|) that \S(h, /i, c)\ < c. 
To understand the case m = 0, one needs deeper results that can be found in Chapter 9 
of |2^] (and were described to us by P. Sarnak). The key result is that the only poles of 
this sum, or of its analogs for other congruence subgroups of SX(2,Z), are at the points 
s + m = 1/2 + it j, where t 2 - — 1 /4 is one of the discrete eigenvalues of the hyperbolic 
Laplacian A = —y 2 (d 2 + d 2 ) on D. A pole at s = 1/2, m = 0, will therefore arise precisely 
if A has a discrete eigenvalue at A = 1/4. It is a happy fact that no such eigenvalue exists 
- the smallest discrete eigenvalue of the Laplacian on D is of order Ai ~ 90. Therefore, 



the sum (|3.52|) may be analytically continued to give a finite value at s = 1/2 for all values 
of m, including the dangerous case m = 0. In addition, simple bounds suffice to show that 
there is no problem with the sum over m. 

This argument shows that all n 7^ Fourier coefficients of the partition function ( |3.1| ) 
are finite. Together with the results of the previous subsection for n = 0, this allows us 
to conclude that the regularization scheme described in Sec. 3.1 provides a finite answer 
for the partition function (|3.1|). We conjecture that, just as for h = 0, the results are not 
compatible with a Hilbert space interpretation. 
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3.5. Aside: the Tree Level Sum over Geometries 

For comparison, we will now consider the sum which arises if we neglect the one-loop 
contribution to the path integral described in Sec. 2.2. As we will see, in this case the 
answer is if anything even worse. 

The sum over geometries, if we take account of only the classical action and not the 
one-loop correction, is 

Z*(t) = ^exp{27r/clm7r}. (3.53) 

c,d 

This sum is quadratically divergent, since at large c and d the summand approaches one. 
Before attempting to regularize the sum, let us compare it to the "correct" sum ( |3.4|) that 



does include the one- loop correction. The "correct" formula has a factor of 1 / VTmr\r]\ 2 
that is outside the summation. This factor is certainly important, but it does not affect 
whether the sum converges. There is also a factor of |1 — q\ 2 inside the sum. This factor 
also turned out in the above analysis to be less important than it may have appeared; we 
simply expanded it as 1 — q — q + qq, and wrote the partition function as a sum of four 
terms. All four terms were similar, and nothing particularly nice happened in adding them 
up. 

The factor in (|3.4|) that actually is important, resulted from the one-loop correction, 



and has no analog in the "naive" sum ( |3.53| ) is the innocent- looking factor of vTmr| 7 that 



is inside the summation. Because of this factor, we had to evaluate our Poincare series at 
s = 1/2. 

This factor is absent in eqn. ( ^.53| ), so now, if we try to define the naive sum Z* 
by introducing a parameter s as before, we will have to evaluate the resulting function at 
s = 0. In fact, the necessary s-dependent function was already introduced in eqn. ( |3.15| ); 
Z*(r) is formally given by the function E(s, k, 0) at s = 0: 

Z*(r) = \imE(s,k,0). (3.54) 

s— >0 

However, as an analytic function in s, E(s, k, 0) has a pole at s = 0. To see this, consider 
the expansion ( |3.34|) of the part of E(s,k,0) which is constant in x. The m = term in 
this sum vanishes, because of the pole in T(s) at s = 0. The m = 1 term gives 

which has a pole at s = coming from the harmonic series £(1) = oo. 
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Hence, without the one-loop correction, the divergence of the sum over geometries 
becomes more serious. Of course, even if we could make sense of the function Z*, we still 
might have trouble giving it a Hilbert space interpretation. 

One might wonder what happens if we modify the definition of Z* to include a finite 
subset of the Brown-Henneaux states - for example, the states of negative energy. This 
means that before summing over geometries, we multiply the exponential of the classical 
action by a polynomial Yin m=o a n,mq n q rn '- The sum over geometries is a sum of terms 
each of which is similar to what was just described, with one such term for each monomial 
q n q m . Each individual monomial will contribute a pole at s = 0, and generically this will 
survive in the sum. 



4. Possible Interpretations 

So far, we have analyzed the sum of known contributions to the partition function of 
pure three-dimensional gravity. As we learned in Sec. 3, the resulting function cannot be 
interpreted as Tr exp(— f3H) for any Hilbert space operator H. 

We will now address the question of how to interpret this result. The most straight- 
forward interpretation is to take the result at face value. Three-dimensional pure gravity 
may not exist as a quantum theory; to get a consistent theory, it may be necessary to 
complete it by adding additional degrees of freedom, and there may be no canonical way 
to do this. 

The other possibility is that some unknown contributions to the partition function 
should be added to the terms that we have evaluated. Here all sorts of speculations are 
possible. We will consider two quite different possibilities. 

4-1- Cosmic Strings 

Known consistent models of 2 + 1-dimensional gravity with negative cosmological 
constant arise from string theory. For example, a famous class of models comes from Type 
IIB superstring theory on AdS3 x S 3 x I, where X is either a K3 surface or a four-torus. 

In these models, the dimensionless ratio k = 1/16G is never a variable parameter, 
but always takes quantized values determined by fluxes that are chosen in the compact- 
ification. (The fact that IjG is not continuously variable is actually |13| a more general 



consequence of the Zamolodchikov c-theorem applied to the boundary CFT.) Moreover, it 
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is always possible to have domain walls across which the fluxes jump. The domain walls 
are constructed from suitably wrapped branes. 

In 2 + 1 dimensions, a domain wall has a 1 + 1-dimensional world-volume and so can 
be viewed as a cosmic string. The existence of these cosmic strings makes the models much 
more unified, as regions with different fluxes can appear as different domains in a single 
spacetime. 

The usual AdS3 x S 3 x X models have supersymmetric moduli, and the values of the 
string tension depend on these moduli. There is a particularly interesting supersymmetric 
value of the string tension at which "long strings" become possible. These are strings that 



can expand to an arbitrarily large size at only a finite cost of energy |p4| , p5| . When long 
strings exist, the energy spectrum is continuous above a certain minimum energy, and the 
partition function Tr exp(—(3H) therefore diverges for all (3. 

The numerical value of the long string tension and of the excitation energy above the 
ground state beyond which the spectrum is continuous are proportional to the jump in 
k = £/16G in crossing the string. The threshold excitation energy is of order 1 (above the 
ground state at energy —2k) if the jump in k is of order 1. 

Since well-established models of three-dimensional gravity have such cosmic strings, 
perhaps they also present in minimal three-dimensional gravity, if it exists. In anything 
that one would want to call pure gravity, the string tension T measured in units of l/£ 2 
must go to infinity as k —>■ oo. Otherwise, the cosmic strings would contribute excitations 
at the AdS3 scale, and one would describe the model as a theory of three-dimensional 
gravity plus matter. We have no idea if the string tension should be proportional to k (as 
one might expect for solitons), to k 1 ! 2 (as for .D-branes), etc. If the jumps in k are correctly 
matched with T, then the strings are long strings and the partition function Tr exp(—f3H) 
that we have been trying to compute in this paper is actually divergent. If the jumps in 
k are smaller than this, then the partition function converges, but to compute it might 
involve corrections from the strings. 

If one thinks that the strings should be long strings, then the requirement T >> £ 2 
means that the jumps in k in crossing strings are much greater than 1. There is some 



tension between this and the proposal in [13] that k can take any integer value. If all 



allowed values of k are connected by strings or domain walls (as in the known string/M- 
theory models), then the fact that the jumps in k are large means that the allowed values 
of k are sparse. 
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4-2. Doubled Sum Over Geometries 

The scenario just described is obviously rather speculative, but at least it has the 
virtue of underscoring our point that the range of conceivable unknown contributions to 
the AdS3 partition function is quite large. 

We will now describe a quite different scenario. To motivate it, we return to the 
formula ( j2.19 ) for the classical action of the basic spacetime Mq ; i: 



I = -Anklmr. (4.1) 

We write more explicitly 

I = 2mk{r - f). (4.2) 

Hence the classical approximation for the contribution of this spacetime to the partition 
function is exp(— i") = exp(— 2ixik{r — f)) = q~ k q~ k . We notice that this is locally the 
product of a holomorphic function of k and an antiholomorphic function, and is globally 
such a product if k is an integer (ensuring that q~ k is single- valued). The one-loop correc- 
tion preserves this factorized form, and therefore the formula ( |2.25|) for the exact partition 
function Zq ; i asociated with Mo,i has the same properties; in fact, Zq,i = Fk(q)Fk(q), with 



F k (q) = q~ k - q")- 1 ■ (4-3) 



n=2 



To the extent that known formulations of three-dimensional gravity are valid, this 
sort of factorization holds for the contribution to the partition function of any classi- 



cal geometry. See [£6|] for a detailed example. The gauge theory description [|T8l , |T9|| of 
three-dimensional gravity gives a natural explanation of this. With negative cosmological 
constant, in Lorentz signature, the gauge group is SL(2,M) x SL(2,M.); the theory is a 
product of two decoupled SL(2,M) theories, associated respectively with left- and right- 
moving modes in the boundary CFT, and this corresponds to holomorphic factorization 
in the Euclidean form of the theory. 



In IJT3I , it was suggested that the exact partition function of pure three-dimensional 
gravity is holomorphically factorized. It was observed that if this is the case, and the 
Brown-Henneaux spectrum is exact until one gets above the classical black hole threshold, 
then the partition function can be determined uniquely. Dual CFT's consistent with the 
necessary spectrum have been called extremal CFT's [p7| . In subsequent work, it has been 
shown that the genus 2 partition function of such a CFT can be uniquely and consistently 
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determined [38], and there has been some work comparing it to what would be expected 



from three dimensions |[29|1 , but on the other hand an interesting but slightly technical 
argument has been given which may show that extremal CFT's do not exist [[JU|. Also, 
it has been argued pH]] that extremal CFT's, if they exist, generally do not have monster 



symmetry for k > 1, in contrast to what happens |32fl for k = 1. 

Now let us discuss holomorphic factorization in view of the sum over geometries. 
Associated to an element 

of SL(2, Z) is a classical spacetime M c ^- Its action is obtained by applying a modular 
transformation to ( |4.2| ): 

J 7 (t) = 27n&; (7T - 7f) . (4.5) 

As usual, 7r = (or + b)/(cr + d), jf = (af + b)/(cf + d). The partition function of the 
manifold M Cj d is 

Z c ,d= F k (q)\ y F h (q)\ y , (4.6) 

and is holomorphically factorized just like Zq^. 

However, when we sum over geometries to evaluate the partition function 

Z= J2 Fk(q)UF k (q)\^, (4.7) 

76W 

holomorphic factorization is lost. (Here W is set of classical geometries M C) d, isomorphic to 
the coset space PSL{2, Z)/Z, with Z being the upper triangular subgroup of SL(2, Z).) In 
fact, this formula is not simply a product of holomorphic and antiholomorphic functions, 
but a sum of such products, somewhat like the partition function of a rational conformal 
field theory. (However, in that case, each term appearing in the sum is separately invariant 
under T : r — > r + 1.) 

Holomorphic factorization has been lost because the sum over topologies is a common 
sum for left- and right-movers. What could be added to restore holomorphic factorization? 
This question has a simple answer, though whether the answer is really relevant to three- 
dimensional gravity remains to be seen. 

If we formally introduce separate topological sums for holomorphic and antiholomor- 
phic variables, defining an extended partition function 

Z= Fk(q)\ 7 F k (q)\ y , (4.8) 

7,7'ew 
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then holomorphic factorization is restored (if the sum converges or can be regularized 
in a satisfactory way). After all, fl4.8|) can be written in the manifestly holomorphically 



factorized form 

K lEW J \ 7 '6W / 

The question (apart from regularizing these sums) is to interpret the double sum over 
topologies. 

What is the classical action corresponding to a given term in the sum in ( f4.9| ) ? We can 
answer this question by simply applying separate modular transformations to holomorphic 
and antiholomorphic variables on the right hand side of eqn. ( f4.2|) . The "action" of a 
"classical solution" that would lead to the term in ( f4.8| ) labeled by a pair 7, 7' must be 

J 7>y = 2iri (jt - 7 ; r) . (4.10) 

For 7 7^ 7', this formula is not real, so it is not the action of a real solution of the Einstein 
equations. 

The most obvious way to try to interpret the formula is to interpret it as the action 
of a complex-valued solution of the Einstein equations. As we noted in the introduction, 
naively speaking the Euclidean path integral is a sum over three-manifolds M that obey 
the boundary conditions. But there is no established way to evaluate the contribution to 
the Euclidean functional integral of a given M unless M admits a classical solution of the 
equations of motion that one can expand around. It is because of this that our first step 
in Sec. 2.1 was to classify the classical solutions; what we classified were the solutions of 
the usual real Einstein equations with Euclidean signature. 

It seems strange to simply ignore three-manifolds that do not admit classical solutions. 
Are their contributions meaningless or zero for some reason? 

One obvious way to slightly generalize the usual framework is to consider the com- 
plexified equations of motion - the Einstein equations R^ u = —Ag^, defined in the usual 
way except with a complex- valued but nondegenerate metric tensor g^ v . If such a complex- 
valued solution is given, one can hope that perturbation theory around it would make sense. 
In the present theory, in which perturbation theory about a classical solution is one-loop 
exact, one could hope to make sense of the contribution of a complex saddle point (that 
is, a complex-valued solution of the equations of motion) to the partition function. 

Unfortunately, we have been unable to find a convincing family of solutions of the 
complexified Einstein equations depending on the pair 7, 7' £ W. Since this is the case, 
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we think another possibility is that the nonperturbative framework of quantum gravity 
really involves a sum not over ordinary geometries in the usual sense, but over some more 
abstract structures that can be defined independently for holomorphic and antiholomorphic 
variables. (A similar idea is expressed in Sec. 3.1 of |[29|| .) Only when the two structures 



coincide can the result be interpreted in terms of a classical geometry. 

"Seeing" The Non- Classical Geometries 

None of this is terribly convincing - no more so, certainly, than the discussion of 
cosmic strings in Sec. 4.1. However, we can describe one fact that offers some slight 
encouragement . 

Let us examine the classical limit of the partition function, by which we mean the limit 
of k — > oo with fixed r. First we will consider the ordinary partition function Z, defined 
by summing over ordinary geometries, and then we will consider the extended partition 
function Z, defined in ( |4.8| ) by a double sum whose meaning is unclear. 

In the classical limit, the sum over geometries is dominated by the geometry that has 
the most negative classical action. In other words, we must pick 7 to minimize the quantity 

J 7 = 2ni('jT - 7f). (4.11) 

We can also write 

I 7 = Re (4ttz7t) = Re (-4ni<yf) . (4.12) 

For generic r, the maximum occurs for a unique 7. 

There actually is a complete democracy between the different choices of 7, because 
they are all obtained from each other by SL(2, Z) transformations. Each 7 dominates the 
partition function (in the large k limit) for r in a suitable region of the upper half plane. 
The thermal AdS manifold M^i dominates in the usual fundamental domain |r| > 1, 
|Rer| < 1/2. Across the arc \t\ = 1, Rer < 1/2, there is a Hawking-Page phase transition 
to a region dominated by M^o, which is the Euclidean black hole. The upper half plane is 
tesselated in phases dominated by different classical solutions. Details are further discussed 
in Sec. 6. 

Now let us consider the extended partition function Z. The action, given in eqn. 
( [1.1(J| ) now depends on the pair 7,7', and is complex. What we want to minimize is the 
real part of the action: 

Re/ 7i7 ' = Re (27ri(^r — j'f)) = Re (2-71^77-) + Re (—2iri-f'f) = Re (271777") + Re (27177V) . 

(4.13) 
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The minimum is always at 7 = 7'. In fact, we can use eqn. (|4 . 1 2| ) to show that 



ReL 



7,7 



-(/ 7 (r)+/ y (r)). 



(4.14) 



Whatever is the minimum with respect to x of T x (r), the minimum of Re/ 7)7 /(r) is at 
7 = 7' = x. 

The conclusion is that, even if non-classical geometries exist and the correct partition 
function is a sum over non-classical geometries as well as classical ones, the semiclassical 
limit of large k is always dominated by a classical geometry. Non-classical geometries exist, 
but they never dominate the semiclassical limit. 

What do we have to do to "see" a non-classical geometry in the semiclassical limit? 
Naturally, we have to ask a non-classical question. The partition function as we have 
defined it so far is Z(r) = Tr exp(27rirLo — 2mfLo), where r is a point in the upper half 
plane IK and f is its complex conjugate. In terms of r = — f , which also takes values in IK, 
we have Z(r) = Tr exp(27rzrLo + 2tivtLq). Now instead of defining r as — f, let us relax 
this condition and simply think of r as a second point in IK.@ Then we should rewrite 
( |4.14|) in the form 



Now it is clear that given any pair 7, 7', we can pick the pair r, r so that the extended 
partition function Z is dominated in the semiclassical limit by 7, 7'. We simply pick r so 
that I x {t) is minimized for x = 7, and r so that I x (j) is minimized for x = 7'. 

Thus, for each non-classical geometry defined by a pair 7, 7', we can find a question 
that that geometry dominates in the semiclassical limit. But to do this, we have to ask a 
rather exotic question that itself depends on a rather unusual analytic continuation. 

5. Black Hole Entropy and its Corrections 

In this and the following section we will discuss a few implications of the assumption 
of holomorphic factorization. In this section we will discuss black hole entropy, and how 
it can be reproduced in holomorphically factor ized theories. 

5 The partition function of a theory in which the high energy density of states is given by 
the entropy of a BTZ black hole remains convergent after this analytic continuation to general 




(4.15) 



t,t e IK. 
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In Sec. 2, we computed the perturbative corrections to the saddle point action for 
the geometry Mq^. Using a modular transformation, this leads to a new formula for 
the subleading corrections to the entropy of the BTZ black hole. We will compare these 
new subleading corrections with those that occur in the extremal holomorphic partition 
functions of fO . 



5.1. Subleading Corrections to the Bekenstein- Hawking Formula 

As we noted in Sec. 2, the geometry we called M^o is just the Euclidean continuation 
of the BTZ black hole. The modular parameter r is related to the Hawking temperature 

and angular potential 6 of the black hole, by r = 6 + if3. The black hole partition 
function is found by applying the modular transformation r — > — 1/r to the expression 
( |2.23|) for the partition function Zq^i of thermal AdS. This gives the one- loop corrected 



partition function for the Euclidean BTZ black hole 

Z lt0 = Z(r)Z(f) (5.1) 
where we have defined the holomorphic piece of the BTZ partition function 

-0-1/24)/., \ 

In this formula we have defined q- = e _27 "/ T . It will also be useful to write Z{r) as 

oo 

Z(r) = C A <£ (5-3) 

A=-fc 

where the coefficients Ca are 

C A =p(A' + k)-p(A' + k-l). (5.4) 

Here p(N) denotes the number of partitions of the integer N. 

Equation ( |5.2| ) encodes various quantum corrections to the thermodynamic properties 
of the BTZ black hole. Equation (|5.2Q is a canonical ensemble partition function, so one 
may, for example, compute the black hole entropy using the usual formula 

S(J3, 9) = log Zi,o - (3 Z^ — ^ (5.5) 
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where r = 9 + i/3. This is a canonical ensemble entropy, evaluated at fixed temperature 
and angular potential 9. 

For the purpose of comparing with CFT predictions, however, it is more useful to 
compute the microcanonical entropy, which counts the number of states N(M, J) at fixed 
energy M and angular momentum J. The energy M and angular momentum J of a state 
are related to left and right-moving dimensions of the state byi 

M = A + A 

x (5-6) 
J = A- A. 

So we may write this density of states as iV(A, A). This density of states is computed 
from the partition function Z\${t, f ) by a pair of Laplace transforms 

/ dr)[ df)q- A q- A Z lfi (r 1 f). (5.7) 

Jie—oo / \Jie — oo / 

These are the usual Laplace transforms that appear when going from canonical to mi- 
crocanonical ensemble. We should emphasize that in this expression r and f should be 
regarded as independent variables. 

We should note that, since ( |5.2|) is a semiclassical partition function around a black 
hole background, it will not have an exact interpretation as a quantum mechanical trace. 
Moreover, the quantization of A and A will not be visible in this approximation^ So 



equation (|5.7| ) should be thought of only as a semi-classical approximation to the number 
of states with dimension (A, A) in the exact quantum theory. Still, we want to explore 
this semiclassical formula. 

The microcanonical entropy S^A, A) is just 

S(A, A) = log N(A, A). (5.8) 

Of course, since Z\$ is holomorphicaly factorized, we may write 

iV(A,A) = A r (A)A r (A), S(A, A) = S(A) + S(A) (5.9) 



6 We are working in AdS units, with i = 1. 

7 To be more precise, in an exact theory with finite entropy, N(A, A) will be a sum of delta 
functions of A, A multiplying positive integers; but the semiclassical approximation does not have 
this form. 
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where 

JV(A) = / drq- A Z(T). (5.10) 



We may use the expansion to get 

00 />ie+oo 

AT(A) = J2 C * / V 

= ^2 Ca drexp I -2ni I Ar H 

A/— i. Jie—oo I V ^ 



(5.11) 



To do this integral, recall the following representation of the Bessel function as a contour 
integral: 

l!(z) = -L iVV z /W)^ (5.12) 

27TZ J 

where the contour encloses the origin in a counterclockwise direction. Taking t 1/t and 
letting 2 = 47ta/— AA 7 , this becomes the integral appearing in (|5.11|) . We end up with the 
following formula for the microcanonical entropy: 



^ / — A' 

AT (A) = e s ( A) = 2n C' A W^-/i(47rv /z AA 7 ) (5.13) 

A'=-k 

where Ca was defined in (|5.4|) . Formula ( p.!3| ) is the main result of this section. 

In the semiclassical approximation, this formula will be dominated by the term with 
A' = —k. We may then use the asymptotic formula for the Bessel function 

h{z) = ~^=e z (l- lz- x + ... ) , atz^oo (5.14) 
\?2nz V ° / 

to get 

13 1 

5(A) = log Af(A) = 4nVkA + -log k - -log A - -log 2 + ... (5.15) 

The first term is the usual Bekenstein-Hawking term, proportional to the area of the BTZ 
black hole. The other terms are logarithmic corrections that typically appear when the 
entropy is computed in microcanonical, as opposed to a canonical, ensemble. They were 



computed for the BTZ black hole in |33[ . 

Of course, equation ( |5.13|) contains many interesting subleading terms in addition to 
these logarithmic terms. Let us consider the terms in (|5.13|) with A > —k. For large values 
of A/k, which is to say for black holes whose radius is large in AdS units, these terms are 
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exponentially subleading. To see this, let us compare the A = —k and A = —k + 1 terms 
in the sum using the asymptotic behavior of the Bessel function ( |5.14p . We find 



Ji(47rV(fe-l)A) 

Thus for A/k large the sum is dominated by the A = —k term plus terms which are 
exponentially small in y/ A/k. For black holes whose size is of order the AdS scale, these 
additional terms can become relevant. 

We should make one more important comment about the structure of the subleading 
terms in (|5.13 ). The terms with A' < are qualitatively of the same form as the leading 



contribution described above. The terms with A' > are qualitatively different: for these 



terms the factors of y— A 7 appearing in ( |5.13| ) are imaginary, and we must worry, for 
example, about which branch of the solution ( |5.13| ) we should take. At this point one 
might wonder about the physical interpretation and implications of these funny terms. 
However, as we will see below, it is precisely these unusual-looking terms which are absent 
in the correct microscopic computation of the black hole entropy. So we will not worry too 
much about their appearance here. 

Our goal is to understand to what extent the microcanonical entropy formula ( |5.13| ) 
can be reproduced in a holomorphically factorized partition function. Of course, we could 
attempt to reproduce the canonical formula (|5.5[) instead, as these two expressions contain 



precisely the same information. But we will find the problem of computing ( |5.13| ) to be 
technically simpler. Our conclusion will be that the terms in (|5.13|) with —k < A' < are 
reproduced in the full partition function. 

Before proceeding, let us ask under what circumstances we expect the formula ( |5.13j ) 
for the entropy to be a good approximation to the number of states of the exact quantum 
theory. Of course, we must take k to be large, in order for the semi-classical approximation 
to be valid. In addition, we must take A to be large so that the corresponding black hole 
is large in Planck units and dominates the entropy of the system. 

5.2. The Rademacher Expansion for Holomorphic Partition Functions 

Let us now consider a modular-invariant, holomorphic CFT with central charge c = 
24k. For such CFTs, k must be an integer, so the holomorphic part of the partition 
function can be expanded as 

oo 
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where the Fa are some positive integers. Modular invariance fixes Z(r) in terms of the 
finite set of coefficients {Fa| — k < A < 0}. Once these polar coefficients are fixed, an 
explicit formula can then be given for Z{r) as a polynomial in the j-invariant, although 
we will not write this formula here. The important point is that all of the coefficients Fa 
with A > are fixed in terms of the coefficients -Fa with A < 0. The formula for the Fa 
with A > in terms of these polar terms is 

F A = 2n £ ^^f^lsiAAAln^f^y^AWl) (5.18) 



where S^A, A',n) is the Kloosterman sum defined in (|3.21|) . This is a convergent series 



expansion of Fa- Expansions of the form ( |5.18| ) are known as Rademacher expansions. 



Rademacher expansions were first applied to three-dimensional gravity in [10 



In [13], a conjecture was made regarding the form of the partition function ( |5.17|) . 



The conjecture is that the theory contains no primary fields with —k < A < 0. CFTs 
with this property are known as extremal CFTs, since in a sense this is the most extreme 
conjecture possible: modular invariance forces one to include additional primaries with 



A > 0. In terms of the coefficients Fa, the conjecture of []T3j is that 



F A = C A for - k < A < 0. (5.19) 



Here Ca are the perturbative coefficients defined in 

The extremal partition function (or any similar function in which the number of 
primary fields of dimension < k is not too large) is consistent with the leading Bekenstein- 
Hawking term in the black hole entropy for A, k large. We are now in a position to ask 
whether subleading terms given by ( |5.13|) are reproduced as well. First, we note that the 
sums ( |5.18| ) and ( p.!3| ) are quite similar. In particular, when n = 1 the Kloosterman sum 
is 1. In this case the summand in ( |5.18j ) coincides precisely with that in ( |5.13| ), at least 
when A' < 0. Thus the Rademacher expansion ( |5.18| ) contains the first k terms in the 
one-loop corrected black hole entropy ( |5.13| ). This is the main result of this section. 

There are two important differences between ( |5.18|) and (|5.13 ). The first is that the 



sum in ( |5.13|) is over — k < A' < oo, while the sum (|5.18|) is only over — k < A' < 0. So 



the funny terms in (|5.13|) with A' > are absent in the full microscopic entropy. It would 
be interesting to understand the physical interpretation of these terms, but we will not 
attempt to do so here. 
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The second important difference is that the Rademacher expansion ( |5.18j ) includes 
terms with n > 1. These terms provide exponentially small corrections to the black hole 
entropy. To see this, we may again use the asymptotic behavior of the Bessel function 
( |5.14|) . For black holes with large A the ratio of the n = 1 and n = 2 terms in the sum 
(|7L8|) goes like 



J -(W A I A 'I> ~ e -W5|*I. (5.20) 

So at large A, the n > 1 terms in the Rademacher expansion are exponentially small in 
A. Again, it would be interesting to study the physical interpretation of these terms. 



6. The Hawking-Page Phase Transition 

In this section we address an apparent puzzle posed by holomorphic factorization and 



the Hawking-Page transition jl5|. As we will see, the resolution to this puzzle is that if 
holomorphic factorization is valid, then the Hawking-Page phase transition is described 



by a condensation of Lee- Yang zeroes of the partition function p4] , |35|1 . We examined the 
zeroes of the partition function as a result of a question from M. Kaneko, and the zeroes we 
find are similar to what have been obtained in investigations |3(],|37| of certain somewhat 
similar modular functions. 

The partition function Z(r) computes a canonical ensemble partition function at fixed 
temperature Im r and angular potential Re r. As described in Sec. 2, for each value of r, 
an infinite family of classical geometries M c ^ will contribute to Z(r). The pair of integers 
(c, d) labels an element of the coset SL(2, Z,)/Z, where Z denotes the shifts r — ► r + n. In 
Sec. 2.2, we computed the contribution of M c ^ to the partition sum Z(r). If we choose 
an element °^ SL(2,Z) corresponding to M Cj d, the contribution to the partition 

function was given by -^o,i( "r+d )- 

For a given value of r, one can ask which geometry M C) d has the largest contribution 
to Z(r) in the semiclassical limit k — ► oo. As explained in ||[10| and in Sec. 4.2, this 
question can be answered by maximizing the classical action — 47rlm (jt) as a function of 
7. If r is in the usual fundamental domain \t\ > 1, |Rer| < 1/2, or any of its translates by 
r — > T + n, then the answer is that the dominant classical solution is Mo,i, which describes 
thermal AdS space. For any r, the dominant classical solution can be found by asking 
which value of 7 G SL{2, Z) will take r into the fundamental domain or one of its translates 
under r — > r + n. For example, M^o, whose Lorentzian continuation describes a black 
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hole in equilibrium with its Hawking radiation, dominates whenever there exists an n such 
that — l/(r + n) lies in the fundamental domain. Roughly speaking, this is a regime of 
low temperature and angular potential. More complicated geometries M c ^, which do not 
have simple Lorentzian interpretations (as they do not admit an analytic continuation to 
a real solution with Lorentz signature) dominate in other regions of the upper half plane 
0~C. Putting this together, we find that the phase diagram of three-dimensional gravity, in 
the semiclassical limit, is given by a sort of subtesselation of the usual tesselation of "K by 
fundamental domains of SX(2,Z). 

We get a subtesselation because certain boundaries between fundamental domains can 
be crossed without any jump in the dominant geometry (an example is that, for Imr >> 1, 
there is no jump in crossing between domains related by r — > r + n). The usual tesselation 
of "K by fundamental domains is sketched in Fig. 3a, and the subtesselation relevant to 
three-dimensional gravity is sketched in Fig. 3b. 



Figure 3: a) The standard tesselation of the upper half r-plane by £X(2, Z) 
fundamental regions, b) The subtesselation that represents the phase diagram 
of three dimensional gravity. The phase boundaries, represented by solid black 
arcs, connect fixed points of <SX(2, Z) of order 3. This characterization was 
found by considering special cases and then using modular invariance. 

It is believed that the phase structure as a function of r of a wide range of three- 
dimensional gravity theories (with different sets of fields) is as we have just described. For 
any finite value of the dimensionless ratio k = £/16G, the partition function is smooth as 
a function of r, but for k — > oo or G — > 0, it becomes non-smooth along the curves drawn 
in Fig. 3b. 

We will address here a puzzle that arises if one assumes holomorphic factorization. 
In this case, we take k to be a positive integer. The partition function factorizes as 
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Z = Zk{q)Zk{q) 1 where the function Z^ is a holomorphic and modular-invariant function. 
Zfc has an expansion near q = involving the negative energy Brown-Henneaux states, 

Z k =q~ k + ..., (6.1) 

and has no other singularities. In order for Z to exhibit the phase structure sketched in 
Fig. 3b in the semiclassical limit, Z^ must do the same. 

The poses a puzzle, since is holomorphic. There is no problem for a sequence of 
smooth but not holomorphic functions depending on a parameter to become non-smooth 
along a phase boundary in some limit. But how can this happen to a sequence of holo- 
morphic functions? 

This question was actually answered by Lee and Yang J51J55 ], The original idea of 



Lee and Yang is that although a system in finite volume can have no phase transition, 
its partition function, in its dependence on the complexified thermodynamic variables, can 
have zeroes. Then, in the infinite volume limit, the zeroes become more numerous and may 
become dense or "condense" along a certain arc, and a true phase transition can emerge. 

In our problem, the limit k — > oo is analogous to a thermodynamic limit, since the 
effective size of AdS space (accessible at a given temperature) grows with k. In a funda- 
mental domain compactified by adding its "cusp," the function Zj. has precisely k poles, 
exhibited in ( |6.1|) . Hence the number of zeroes in a fundamental domain is also equal to 
k. Here in the case of zeroes that are on the boundary of the fundamental domain, one 
must count only half the zeroes (or one half of the boundary) to avoid double-counting. 

Before going into details, we will give an informal explanation in which we ignore the 
one-loop correction. We consider the phase boundary along the arc \t\ = 1, |Rer| < 1/2 
separating the thermal AdS phase from the Euclidean black hole. We will denote the 
portion of this arc with Re r > by C: 

C = {t = x + iy\\r\ = 1,0 <x < 1/2}. (6.2) 

The other portion of the arc, with Re r < 0, is the image of C under the modular 
transformation r — > — 1/r. Now, the contributions of Mq ; i and Mi o to the holomorphic 
Poincare series (ignoring the one-loop correction) are exp(— 2-Kikr) and exp(— 2nik(— 1/r)). 
On the arc C, we have 1/r = f and the sum of these two contributions is 

exp(— 27rikr) + exp(2nikf) . (6-3) 
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The sum vanishes if exp(2nik(r + f)) = —1. Setting r = exp(z#), n/3 < 9 < 2n/3, we 
need exp(47ri/e cos 9) = — 1 or 47r/ccos# = (2n + l)7r, n G Z. On the arc C, we have 
< cos# < 1/2. So this equation has /c roots. In the limit k — > oo, these roots become 
dense everywhere along the arc. This is the condensation of Lee- Yang zeroes of the partition 
function. 

Now we give a more precise account including the one-loop correction. However, the 
result does not depend on all the details of the function Z^. We will show, following |]36|,|37 



that for a certain class of holomorphic functions, all zeroes lie on the arc C or one of its 
images under SX(2,Z), as shown in Fig. 3. In the large k limit, the zeroes become dense 
on C. 

We consider a modular-invariant function that is regular except for a pole of order 
k at q = 0. Such a function has a Laurent expansion 

oo 

Z k(r) = Fa( ? A - (6-4) 

A=-k 

We require that = 1 and that the coefficients Fa, A < 0, do not increase too quickly. 
More precisely, we will show that if 

F A < e 2-(o.6i)(A +fc ) 5 for A < (6.5) 

then the zeroes are all as just described. The origin of the numerical coefficient will be 
described below. This bound has a simple physical explanation: when there are too many 
primary fields of low dimension, i.e. too much light matter in the theory, the formation of 
stable black holes is not possible. Problems with black hole physics when there are many 
light species of matter have been considered in, for example, | ]38|| . 

We should emphasize that the extremal partition functions proposed in [[13[ obey 



the bound ( |6.5|) . So in the large k limit the extremal partition functions exhibit the 
condensation of zeroes described above. In fact, for A < and large k the coefficients of 
the extremal partition function grow much more slowly than ( |6.4|) . For large A + k, they 
grow like 

F A ~ e f <K e 2 7r( o.6i)(A+ fc )_ (6 6) 

So one could add many more additional primaries with dimension A < without spoiling 
the phase transition described above. 

Our proof of the foregoing assertions relies on a few specific facts about modular 
functions, which we now review. 
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6.1. Review: Properties of J 

First, we recall a few properties of the modular invariant J(r). Proofs of many of 
our assertions in this and the next section can be found in ||39|| . The J function has a 
q— expansion 

J(r) = V c(m)q m = - + 196884g + . . . (6.7) 
^— ' q 

m>-l * 

where the coefficients c(m) are positive integers. Our choices of overall normalization and 
constant term c(0) = agree with [|32|] but differ from the most classical convention. As 
the c(m) are positive, it follows that J obeys 

J(r) = J(-f ) (6.8) 

and so is real along the imaginary r axis. In particular, it can be shown that J decreases 
monotonically along the imaginary axis from J(ioo) = oo to a minimum at r = i with 
J(i) = 984. (Continuing past i on the imaginary axis, J grows again in view of its 
invar iance under r — > — 1/r.) 
The modular domain 

D = {t = x + iy\\r\ > 1,-1/2 < x < 1/2} (6.9) 

is a fundamental region for SX(2, Z). On the interior of D, J takes every complex value 
exactly once. (This follows from the fact that in the compactified fundamental domain, 
J has precisely one pole, which is at r = zoo.) Using the invariance of J under SX(2,Z) 
and ( |6.8|) , one can show that J(r) takes real values on the boundary dD. In particular, 
it turns out that on the arc r = e l6 *, J(r) decreases monotonically from 984 to —744 as 
9 runs from 7r/2 to n/3. This is the arc we called C above. On the arcs x = ±1/2, J 
decreases monotonically from —744 to — oo as y runs from \/3/2 to oo. 

Although the exact form of J(q) is quite complicated, in many cases the "tree-level 
approximation" 

J~g _1 + ... (6.10) 

is very useful. To this end, we note that for any given value of r = x + iy, we have the 
bound 

\J(t) -q- 1 \ = \J2 c{m)q m \ < ^ c(m)\q\ m = J(iy) - e 2?ry . (6.11) 

m>l m>l 
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The function J2 m>1 c(m) \q\ 171 depends only on y and is a monotonically decreasing function 
of y, so in the fundamental domain D, it is bounded above by its value at y = a/3/2, which 
is the minimum of y in D. This value is 

M = J(za/3/2) — e 71 ^ ~ 1335. (6.12) 

So in D, we have 

\J(r) -q' 1 ] < M. (6.13) 

Applying the triangle inequality, it follows that, throughout D, we can bound the value of 
J by 

\J(r)\ < e 2iry + M. (6.14) 

The approximation ( |6.10| ) is a good one when the imaginary part of r is large, i.e. 
when r is close to the cusp at r — ► oo. By considering modular transformations of the 
equation ( |6.13|) we arrive at other approximations which are good in other regions of the 
upper half plane. In particular, for any element 7 G SL(2, Z), if 77- is in the fundamental 
domain D we have 

|J(r) - e - 27 " 7T | < M, (6.15) 

and therefore 

\J(r)\ < e 2wIm7T + M. (6.16) 

This final bound is particularly useful, for the following reason. For a given value of r, 
consider the set of possible values of Im 77", for all 7 G SX(2,Z). This set has a finite 
maximum value which occurs when 7T is in the fundamental domain or one of its images 
under the map r — > r + n. So, using the formula ( |3.9| ) for Im 77", we conclude that for any 
value of r = x + iy in the upper half plane we have the bound 

, / \ 1 f 27r Im r 1 . 
I J(r)| < exp <^ max C)dl — ^ ^ + M. (6.17) 



|cr + d\ 

Here the maximum is taken over all relatively prime c and d. 
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6.2. Properties of T n J 

We will now consider the action of the Hecke operator T n on J, defined by 

5-1 

T - J = E E J (( nr + ^)A 2 )- ( 6 - 18 ) 

5|n (3=0 

This is a new modular invariant function with an n-th order pole at q = oo. Our normal- 
ization (which differs by a factor of n from that in the literature) has been chosen so that 
the q— expansion of T n J starts with q~ n : 

oo 

T n J= J2 c n (m)q m = q- n + 0(q). (6.19) 

m= — n 

The Fourier coefficients c n (m) of T n J are related to the coefficients c(m) of J by 

c n (m)= ^ ^c(mn/5 2 )q m . (6.20) 

8\(m,n) 

From this it follows that the Fourier coefficients of T n J are positive integers, just like those 
of J. 

From these facts, it is straightforward to see that T n J satisfies many of the same prop- 
erties as J. In particular, T n J is real along the imaginary axis, as well as on the boundary 
of the fundamental domain dD. Moreover, just as with J, the "tree-level approximation" 

T n J~q~ n + ... (6.21) 

is good enough for many purposes. This approximation is very good when r — > ioo, but 
actually it can be extended to understand the zeroes of T n J, as was shown in P7| . We 
will summarize their argument here. In section 6.3, we apply the results derived below to 
more general modular-invariant partition functions. 

The basic idea is to apply either ( |6.15| ) or (|6.17| ) to each term in the sum (|6.18| ) . This 



allows us to bound the value of T n J(a) for a point a = x + iy on the arc C defined by 



|er| = 1 and < x < 1/2. Let us first consider the 5 = 1 term in (|6.18|) , which is J(na). 
For any a in C, we can find an integer m such that na + m lies in the fundamental domain. 
This allows us to apply (|6.15|) with r = no and jt = no + m to get 



\J(no) — e ~ 27 ™ ncr | < M. (6.22) 
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Now, consider the 5 = n, (3 = term in ( |6.18| ), which is J(a/n). Since a lies in C, we can 
find an integer m such that —n/a + m lies in the fundamental domain. So we may apply 
( |6.15|) , with r = a/n and 7 = —n/a + m to get 

|J(a/n) -e- 2 " m/,T | < M. (6.23) 

Using these two equations, we may apply the triangle inequality to (|6.18|) to get 



t n J{a) — e - e 1 



(6.24) 



Here the sum ' is over the same set of (3 and 5 as in ( |6.1(j| ) , except we have dropped the 
two terms with 5 = 1 and 8 = n, (3 = 0. We will now apply (|6.17| ) to each term on the 
right hand side of ( |6.24j ). Consider first the term with S = n, (3 = n — 1, which is J f 2 ^). 



We would like to apply (|6.17| ) with 



T 



(7-1 



so we must ask what possible values |cr + d\ 



can take for this value of r. Now, since a lies on C, it follows that \a — 1| > 1 and hence 
|t| = l^-^-H > l/n. This implies that |cr + d\ > 1/n for all possible choices of c and d. 
Since Im r = y/n, equation (|6.17|) gives 

a + n — 1 



J 



n 



< e 2irny + M. 



(6.25) 



Let us consider the case where 5 = n and < j3 < n — 1, where we may apply ( |6.17| ) with 
r = -2-±£. For this range of 0, \a + (3\> y/2, so |t| > y2/n. Hence |cr + d\ > \p2jn and 
( |6.15|) gives 



J 



a + /3 



< e nny + M. 



(6.26) 



Finally, we consider the cases where 1 < 5 < n. In this case we apply ( |6.17| ) with r = 
(na + (36) /S 2 . The fact that a lies on C implies that \cr + d\ > \^3n/5 2 . So we end up 
with the bound 



J ( — ) <e — y +M. (6.27) 
Putting this all together, equation ( |6.24| ) becomes 

|T n J(a) - e" 2 ™ - g-Svrm/a | < + ^ny + rf^y + n ^ M _ ( 6 28 ) 

The factors of n and n 2 in ( |6.28| ) come from the simple fact that there are less than n terms 
with 6 = 71, and less than n 2 terms with 1 < 6 < n. Multiplying both sides of ( |6.28| ) by 
e _27rny , and using the fact that e - 27rm / <T = ^ina f or p j n t g on th e curve C, this becomes 



\T n Je- 2nny - 2 cos(27rnx)| < 1 + ne~™ y + n 2 e~^ y + Mn 2 e~ 2 ™ y . 



(6.29) 
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For the moment, we concentrate on the case n > 2. Since y > \/3/2 for any point on C, 
we may evaluate this right hand side of ( |6.29|) to get 

\T n Je- 2nny - 2cos(27rnx)| < 1.12. (6.30) 

This formula is valid for any point on the arc C and n > 2. 

Equation (|6.30| ) is an important result. It places strong constraints on the location 
of the zeros of T n J. To see this, note that on each of the n intervals I m defined by 
— < x < , m = 0, . . . , n — 1, the function 2 cos(27mx) varies monotonically from —2 
to 2. The right hand side of ( |6.30| ), however, is less than 2. This implies that T n J must 
have n distinct zeroes on C, with one in each interval I m . As n —>■ oo, the lengths of these 
intervals vanish and the zeroes of T n J become dense on C. Finally, note that, as T n J is an 
n th order polynomial in J, it has only n zeroes on D. In particular, it has no additional 
zeroes beyond those described above. 

We have only proven the bound ( |6.30|) for n > 2. For our application, we also need the 
n = and n = 1 cases, which we will consider separately. For n = 0, we define TqJ = 1. 
So the bound ( |6.30|) is trivial. For n = 1, we have T\J = J. In this case we have the 
slightly weaker bound 

\J(r)e- 27Ty - 2cos(27rx)| < 1.22 (6.31) 

for points on C. This is straightforward to verify numerically, although it may also be 
understood analytically. The point is that the function J(r)e _27ry — 2 cos(27rx) is monotonic 
along the arc C. So the value of this function on C is bounded by its values at the endpoints 
r = % and r = l/2 + y/3/2i. Using J(i) = 984 and J(l/2 + y/S/2i) = -744 we may compute 
the values of the function at these endpoints, giving the bound Q). 

6. 3. Zeroes of Modular Invariant Partition Functions 

We are now ready for the general case of a modular invariant partition function at 
level k: 

oo 

Z k (r)= F ^ A = E F aT| A | J(r) (6.32) 

A=-fc A=-k 

where the Fa are non-negative integers and F_k = 1. As with J(r), the fact that Z^{t) 
has real coefficients and is modular-invariant means that it is real on the imaginary r axis 
as well as on the boundary of the fundamental region dD. 
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In the previous section, we proved the estimate ( |6.30|) for T n J on C with n > 2. Along 



with the estimate ( |6.31| ) on C, this implies that for all n > 

T n J = e 2lzny (2 cos(2nnx) + E n ) (6.33) 
where E n is an error term obeying 

\E n \ < 1.22. (6.34) 

By simply adding up the inequalities (|6.33| ) for n = 0, . . . , k, with coefficients F- n , and 
using the fact that F_ k = 1, we get a a similar estimate for Z k : 

o 

Z k e~ 2nky - 2 cos(2ta) = E k + ^ F A e- 27T{k+A)y (2 cos(2nAx) + E ]A] ) . (6.35) 

A=-fe+l 

To bound the location of the zeroes of Z k , we must show that the right hand side is 
less than 2. Then, as in the previous section, the zeroes of Z k will lie on C and become 
dense in the large k limit. 

For this to be the case, F A must not grow too quickly with A. For example, assume 

that 

Fa < Ae 27Ta{A+k) for - k < A < (6.36) 
where a and A are positive constants. In this case 



^2 F A e~ 27T{k+A)y < A ^2 e 27r{k+A)(a ~ y) 

A=-k+l A=-k+l 

k 

< ^ ^ e -27rn(a-V3/2) < ^_ 



(6.37) 



e 27r(a-V3/2) _ 1 ' 
n=l 

In the second line, we have used the fact that y > V3/2 on C. Since 

\2 cos(27rAa;) + E\ A \ \<2 + \E\ A \ \ < 3.22, (6.38) 
it follows from (|6.34j ) that 

\Z k e- 2nky -2cos(2knx) \ < \E k \+3.22 f , (6.39) 

g27r(vo/2 — ct) ^ 
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which is less than 2 for certain values of A and a. In particular, using \E^\ < 1.22 and 
setting A = 1, we find that the right hand side is less than 2 provided that 

a < 0.61 (6.40) 

We have used several rather conservative estimates, so it is quite possible that the zeroes 
of Zk condense on C even when the coefficients .Fa do not satisfy the precise bound given 
above. 

It is easy to check that this bound is indeed satisfied for the extremal partition func- 



tions of [13]. In fact, the coefficients Z&(A) grow much more slowly than the required 
behavior, indicating that many more primary fields could be added without spoiling the 
existence and nature of the phase transition. 



7. Supergravity Partition Functions 

Our goal here will be to repeat the analysis of Sees. 2 and 3 for supergravity. We 
will consider only the basic case of M = 1 supergravity, in which the symmetry group 
SL(2,R) x SL(2,R) of AdS 3 is replaced by OSp(l\2) x OSp(l\2), where OSp(l\2) is a 
supergroup whose bosonic part is Sp(2, M) = SL(2,WL). The boundary CFT has (1,1) 
supersymmetry, that is, M = 1 supersymmetry for both left- and right-movers. 

7.1. Formalism 

In supergravity, there are a few closely related choices of possible partition function. 
We can compute either Tr exp(—(3H — i6J) or Tr (— 1) F exp(— (3H — i6J). And we can 
compute this trace in either the Neveu-Schwarz (NS) or the Ramond (R) sector. As usual, 
we attempt to compute these partition functions by summing over three-manifolds M 
that are locally AdS3 and whose conformal boundary £ is a Riemann surface of genus 1. 
The four possible partition functions (NS or R, with or without an insertion of (— 1) F ) 
correspond to the four spin structures^ on E. The three that correspond to odd spin 
structures are permuted by the action of SL(2, Z), as we further discuss below. 

8 In fact, using a discrete i?-symmetry that appears to be present in the boundary CFT, we 
could choose different spin structures for left- and right-movers on X, but we will not consider 
this generalization. 
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Once we pick the spin structure on E, we then sum over choices of M such that 
the given spin structure on E does extend over M. For example, what spin structure on 
E is compatible with taking M = M 0) i to be the three-manifold related to perturbative 
excitations of AdS3? 

In M 0j i, the "spatial" circle on E is contractible. This means that the NS spin 
structure on the spatial circle extends over M 0) i and the R spin structure does not. Hence 
Mo,i contributes to traces in the NS sector, not the R sector. 

What is the spectrum of thermal excitations in the NS sector? In ordinary gravity, 
the thermal excitations of left-movers are obtained by acting on the ground state |fi) with 
Virasoro generators L_ n , n > 2. When the boundary theory has M = 1 supersymmetry, 
we can also act on \Q) with superconformal generators G_ n+ i/2, n > 2. (We recall that 
GLi/ 2 |0) = 0.) Writing —k*/2 for the ground state energy, the partition function of 
left-moving excitations is therefore 

00 1 _l n n-l/2 

Including both left- and right-moving excitations, the contribution of M 0j i to F(q,q) = 
Tr NS exp(-/3H - id J) is 

2 

(7.2) 

This formula can be justified exactly as in Sec. 2.2 for the bosonic case. 

The complete function F, or more precisely the sum of known contributions to it, 
is evaluated by summing over all those modular images of M^i over which the relevant 
spin structure extends. We can represent the four spin structures on the two-torus E by a 
column vector 

{:)■ p- s > 

where \i and v represent respectively the fermion boundary conditions in the "time" and 
"space" directions on the two-torus E. We consider jj, and v to be valued in ^Z/Z, taking 
the values 1/2 for antiperiodic (NS) boundary conditions and for periodic (R) ones. An 
element of SL(2, Z) acts on the spin structures by 




^b,i — 



-k*/2 



n 

n=2 



1+ n-1/2 



1 - q r 
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The NS partition function F(r) = Ttns exp(—(3H — i9J) corresponds to fx = v = 1/2. 
This condition is invariant under the subgroup of SL(2, Z) characterized by saying that 
c + d and a + b are both odd. If c + d is odd, we can make a + b odd by adding to (a, b) 
a multiple of (c, d). F, or at least the sum of known contributions to it, can therefore be 
computed by summing F 0j i over modular images with c + d odd: 

F(r)= ^d(r) (7.5) 

c,d| c+d odd 



or equivalent ly 



F(t)= F 0tl ((aT + b)/(cr + d)). (7.6) 

c,d| c+d odd 



For any given c, d such that c + d is odd, we pick a and b such that ad — bc = 1 and a + b is 
odd. Because Fo,i(r) is invariant under r — ► r + 2, the summand in ( |7.6[ ) does not depend 
on this choice. 

It is also of interest to compute partition functions with other spin structures. How- 
ever, the three even spin structures on E (the ones for which \x and v are not both zero) 
are permuted by SL(2,Z) and so the associated partition functions are not really inde- 
pendent functions. If we set \i = 0, v = 1/2, we get G(r) = Ttns (—l) F exp(—(3H — i6J). 
The contribution of M^i to this partition function is obtained by reversing the sign of all 
fermionic contributions in ([7.21): 



Go,i(r) 



q 



i n n-l/2 2 



n i _ q > 

n=2 * 



(7.7) 



The subgroup of SL(2, Z) that preserves the conditions = 0, ^ = 1/2 is characterized by 
requiring that 6 should be even, which implies that a and d are odd. Hence 

G(r)= ^ G Cjd = ^ Go,i((ar + 6)/(cr + d)), (7.8) 

c,d|dodd c,d|dodd 

where for given c, (i, we pick a, 6 so that ad— be = 1 and 6 is even. A modular transformation 
T : r — > r + 1 exchanges (fi, v) = (0, 1/2) with (//, z/) = (1/2, 1/2), so in particular 

F(t) = G(t + 1) = F(t + 2), (7.9) 



and the summand in (|7.8| ) does not depend on the choice of a, b. 

60 



The Ramond partition function K = TYr exp(— f3H — i6J) is computed from \i = 1/2, 
v = 0, so 

K(t) = G(-1/t). (7.10) 

This completes our characterization of three of the four partition functions. Finally, in 
any supersymmetric theory with discrete spectrum, the fourth partition function / = 
Ttr( — 1) F exp(—(3H — i9J) is an integer, independent of (3 and 9, since it can be interpreted 
as the index of a supersymmetry generator. It must be computed using the odd spin 
structure, the one with \i = v = 0. In three-dimensional gravity, assuming that the 
partition function can be computed by summing over smooth three- geometries, / vanishes, 
since the odd spin structure does not extend over any three-manifold with boundary S. 

General Structure 

Some simple but important facts follow from (|7.9|). Since F(r) = Ttns exp(—(3H — 
iOJ) is invariant under r — > r + 2, it follows that all eigenvalues of J = Lq — Lq take 
values in Z/2. The transformation r — > r + 1 acts as 1 or —1 on states with integer or 
half-integer J; thus it acts as (— 1) 2J . On the other hand, if we insert a factor of (— 1) F 
in the trace, we get G(r) = Tr NS (-1) F exp(-0H - i9J). Since G(r) = F(t + 1), the 
operator (— 1) F is equivalent to (— 1) 2J . In other words, states of integer or half-integer 
J are bosonic or fermionic, respectively. The exact theory (to the extent that it can be 
reconstructed from the sum over smooth classical geometries) inherits this property from 
the perturbative spectrum of Brown-Henneaux excitations. 

The general form of G(r), assuming that it has a Hilbert space interpretation, should 
therefore be 

G(r) = ^2^2 a n ,j exp(-PE n> j + i9j) - XI S ^ ex P(-/^n,j + ie 3)- (7-11) 

jei n jeZ+1/2 n 

Here E n j, n = 1, 2, 3, . . ., are the energy eigenvalues for states of angular momentum j, 
and a n j is the number of states of energy E n j. The a n j are positive integers; the minus 
sign preceding the second term in ( [7.1 1| ) reflects the relation (— 1) 2J = (— 1) F . 
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7.2. The Computation 



We we will now compute the partition functions of the previous section by summing 
over smooth geometries. We will consider the partition function G(r) = Ttns ( — 1) F exp(—(3H- 
iOJ), as it is technically the simplest; the two other non-zero partition functions are then 
given by O and (fMOft . 

From (|7.7|) and (fT.81) we have 



c,d| d odd 



n 

n>2 



1 _ g»-l/2 



1 — (f 



(7.12) 



To understand the modular transformation properties of this sum, it is useful to rewrite 
the infinite product in terms of Dedekind eta functions. Using the identities 



11(1-5' 



71 = 2 



and 



this may be written as 



n(i-? n - 1/2 : 



n=2 



.! _ q^jl - q) 
t]{t) 

_ q 1 / 48 t]{t/2) 
"(l-gVa) rj(r) 



(7.13) 



(7.14) 



g(t) = y ( ir fc * +3 / 24 ii + q 1/2 \ 2 



(7.15) 



We may now extract these Dedekind eta functions from the sum, using the fact that 



Vim r 1 77 (t ) 1 2 is modular invariant: 

Wt/2)\ 2 



G(t) 



E(?/ 1/2 kr fc * +3/24 ii-9 1/2 

(E(k* - 3/24, 0) + E(k* + 1 - 3/24, 0) + 



y 1/2 \v(r)\ 4 

\v(r/2)f 

y 1/2 \v(r)\ 4 

E(k* + 1/2 - 3/24, 1/2) + E(k* + 1/2 - 3/24, -1/2)) 



(7.16) 



In the second line we have defined 



E(k,h)= 



y 



1/2 



■ exp {2itk Im 7T + 2m n Re 7r} . 



c,d\ d odd 



|ct + d\ 

This is the supersymmetric version of the Poincare series studied in Sec. 3. 



(7.17) 
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This sum is divergent, for the same reasons described in Sec. 3.1. In particular, at 
large c and d the exponential approaches one and we are left with the linearly divergent 
sum Y2 C d |cr + (i| _1 . The sum may be regularized by considering the more general Poincare 
series 

£'(s,k, / u)= "j — -j— exp {27tk Im 7r + 27ri// Re 7r} (7-18) 

c,d\ d odd 

as an analytic function of s and taking s — > 1/2. This regularization scheme can be justified 
on physical grounds, following the same line of argument presented in Sec. 3.1. 

We will now proceed to compute the Fourier coefficients of the sum ( [7.18|) . Our 
calculation is quite similar to that done in Sec. 3. The only differences are that in ( |7.18|) , 
k and \x are allowed to be half-integer, and that we are summing over c and d such that 
(c, d) = 1 and d is odd. These two conditions on c and d can be combined in to the single 
condition (2c, d) = 1. 



We start by letting d = d' + 2nc, where d' = d mod 2c. The sum in (|7.18|) can be 
written as a sum over c, d! , and n: 



E(s 1 K ll ,) = y s e 2 ^y+^ + J2 E E/( c ' rf » ( 7 - 19 ) 

c>0 d'£Z/2cZ n£Z 



where 



, , 2,71 Ky fa cx + d 

f(c, d,n)= — exp + 2ttz^ - 



|c(t + 2n) +d'\ 2s { |c(r + 2n) +d'\ 2 c\c{r + 2n) + d'\ 2 

(7.20) 

We will now apply the Poisson summation formula to the sum over n, as we did for the 
bosonic partition function in Sec. 3.2. First, we must compute the Fourier transform 



/oo 
dn e 2 ™f{c,d',n) 
-oo 



= I exp (tojW> a -™ - nin X ) f dt (* \ * exp ' 2 ^ y ~ ^ 



2 2c 7 J-oo Vc 2 (t 2 + y 2 )7 I c 2 (t 2 + y 2 

(7.21) 

We have written the integral in terms of a shifted integration variable t = 2n + x + — . 



Recall that, as described after equation (777), we must chose our integers a and b such 



that b is even. The determinant condition ad— be = 1 therefore implies that ad' = 1 mod 2c. 
Since 2\x G Z, this implies that the exponential prefactor in (|7.21| ) is a function only of 
d' = d mod 2c. We may therefore extract the d! dependence of the sum ( |7.18|) into the 
Kloosterman sum S{— n, 2/i; 2c), as defined in ( |3.21|) . The integral appearing in ( |7. 2 1|) is 
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precisely that defined in ( |3.24j) . The Fourier cofficients of the Poincare series (|7.18 ) are 
therefore given by 

E(s, k, fi) = y s e 27T ^ y+l ^ + e™^(s, k, fi) (7.22) 

n 

where 

-. oo / oo \ 

E^s, «, M) = 2 E / m > n/2( S > K ' /*) y 1_m ~' E c" 2(m+s) S(-n, 2 M ; 2c) (7.23) 

m=0 \c=l / 



is defined in terms of the integrals ( |3.24| ) . 



We will now restrict our attention to the n = case. First consider [i = 0. In this case, 



the integrals were given in (|3.26|) . To do the sum over c, note first that the Kloosterman 



sum S(0, 0, 2c) is equal to Euler's totient function </>(2c). The sum over c is 

c>0 c>0 SV ^ ~ J) 

To prove this formula, we first evaluate the sum J2 C odd C_<T( ^( C )' wnere we have defined 
a = 2(m + s). Multiplying this sum by 

E = (1 - 2-)C(<r) ( 7 - 25 ) 

n odd 

gives 

(i - 2-')c(«t) x; c_cr ^ c ) = e m-^w 

c odd c,n odd 

= E m_,J (e^ c )) • ( 7 - 26 ) 

m odd \ c|m / 

= £ m-+ 1 = (l-2 1 -)C(a-l) 

m odd 

In the second line we defined the new variable m = cn to be summed over, and in the third 
line we have used the basic identity for the totient function J2 c \m ^( c ) = m - We thus have 

V c-m = f" 2 ^" 1 '. (7.27) 

The sum over even c may be done by recalling the identity that was proved in Sec. 3 by 
similar methods: 

£ c - W= «£L_i2 (7 . 28) 

C>1 ^ ' 
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to get 

T c-^(c) = 1 C{tT ~ 1] . (7.29) 



c even 

But 



c~^(c) = 2-*£y*#2c). (7.30) 

c even c>0 

Combining the last two formulas, we arrive at ( 7.24]) . 

As in the bosonic case, we must be careful when taking s — > 1/2. For m = (and 
n = 0), the sum ( |7.24p vanishes at s = 1/2, whereas the integral Jo,o has a pole at 
m = 0, s = 1/2, as we see in ( |3.26| ). The product of the two factors is finite; in fact, 
r(s — l/2)/C(2s) — > 2 at s — > 1/2. The m > terms are finite without any such subtleties. 
Evaluating the first three terms in the expansion of ( |7.18| ) gives 

g o(1 /2, K , 0) = -y^ + ( J*L„) „- V + (jf^) + (7.31) 

Let us now consider the /i = ±1/2 terms. In this case the integrals are only slightly 
different from those done in Sec. 3, which were evaluated at \i = ±1. For m = 0, we find 
that 

J 0| o(a, «, ±1/2) = v^ r(g r ~^ /2) ■ (7-32) 

This is the analog of equation (ft.28Q used in the computation of the bosonic partition 
function. For m > the integrals are complicated hypergeometric functions of the sort 
written down in (|3.27| ). At s = 1/2 these hypergeometric functions simply to give 

21-77-1^771+1/2 

/ m ,o(l/2, k, ±1/2) = T m (2K). (7.33) 

ml (m + 1/2) 

This is the analog of equation ( |3.42| ) used in the bosonic case. 

To do the sum over c, we use the fact that the Kloosterman sum S(0, ±1, 2c) is equal 
to the Mobius function //(2c). The sum over c is given by 

n2(m+s) 

y c- 2( - m+a ^S{0, ±1, 2c) = V c" 2 ( m + s V(2c) = - - 9 , - , — — -. (7.34) 

Z_. ^ ' ^ > (2 2 (-+ s )-l)C(2(m± S )) 1 J 

We can obtain this formula as follows. We first compute the sum J2 C odd c _ °"/i(c). Multi- 
plying this sum by ( |7.25| ) we get 



odd 




(7.35) 
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In the second line we have let m = cn and in the third line we have used the fact that 
for any m, X] c |mM c ) = &m,i- The sum over even c may be done by recalling the similar 
identity from Sec. 3 

E C "W = 7^T (7 ' 36) 

01 ^ ' 



to get 



£ «=-««=) = ~ikim <7 - 37) 



c even 



which is equivalent to (|7.34|) . 

At s = 1/2, we again must be careful to cancel the zero in ( |7.34j ) at m = against a 
pole in (|7.32| ), using the fact that T(s — l/2)/£(2s) — > 2 as s — > 1/2. Including the next 
two terms in the series, we find 

E (l/2, k, ±1/2) = -2,V2 _ ,-1/2 _ (^(8 K 2 - 1)) + 0(y~^) 

(7.38) 

Putting this all together gives the following expansion for the partition function 

GW ^' W + W^ 6+ 21«3) * 

47r 2 (2880A;* 2 -16A;*(27r 4 - 135)+45 - 127r 4 ) _ 2 ^. 1 ' ' 
4185C^ V ~ + ° iy ~ ] J- 

As in the bosonic case, a Hilbert space interpretation is precluded both by the fact that 
the coefficient of the leading correction to Gro,i( r ) is negative and by the fact that there 
are additional corrections involving powers of 1/y. A Hilbert space interpretation would 
require the leading coefficient to be positive in view of ( |7.11| ) ; the leading correction to the 
Brown-Henneaux spectrum arises for an integer value of the angular momentum, namely 
j = n = 0. 

We have so far considered just the n = case. The discussion for the n ^ case will 
be qualitatively similar to the bosonic case discussed in Sec. 3.4. The integrals can be 
evaluated explicitly - they are the same ones we considered in Sec. 3.4 - and contribute 
terms to the NS sector partition function which fall off exponentially at large y. As in the 
bosonic case, the sums over c cannot be expressed in terms of elementary number theoretic 
functions but are related to an appropriately defined Selberg zeta function. In particular, 
they should be related to the Selberg zeta function associated to the congruence subgroup 
described in Sec. 7.1, where b is constrained to be even. As in Sec. 3.4, this should provide 



66 



a finite regularization of the sum because the associated Laplacian does not have a discrete 
eigenfunction with eigenvalue 1/4. 

We conclude this section by emphasizing a few possible interpretations of the result 
( [7.39| ) for the supergravity partition function G(r). As in the bosonic case, this partition 
function does not have the structure of a Hilbert space trace Tr NS (-l) F e -^ H - l9J . One 
possible interpretation of this result is that three-dimensional pure supergravity does not 
exist as a quantum theory, and that additional degrees of freedom must be included in order 
to render the theory sensible. A second possible implication is that the pure supergravity 
theory exists, but contains long strings, similar to those described in Sec. 4.1. A third 
possibility is that additional complex saddle points must be included in the sum over 
geometries, as described for the bosonic theory in Sec. 4.2. This might lead to a sum 
over two copies of the modular group, resulting in a holomorphic partition function which 



might coincide with the extremal supergravity partition functions described in [|T3 



7.3. Phase Transitions 

In this section we will comment briefly on the supergravity generalization of the results 
of Sec. 6. In particular, we ask whether, for holomorphically factorized supergravity 
partition functions, Hawking-Page phase transitions occur via a condensation of zeroes 
along curves in the complex temperature plane. Numerical evidence described to us by M. 
Kaneko suggests that this may be the case. 

Let us first ask what the phase diagram of three-dimensional supergravity should be. 
We will consider the NS partition function F(r) = Tr^s e~ l3H ~ l9J . As described in Sec. 
7.1, F(r) is invariant under the subgroup of SL(2, Z) defined by the condition that a + b 
and c + d are both odd; this group is sometimes called IV The tessellation of the upper 
half plane by fundamental domains of Tg is shown in Fig. 4a) . 





Figure 4: a) The tessellation of the upper half r-plane by fundamental 
domains of Tg. b) The subtessellation of the upper half plane corresponding 
to the coset Tg /Z (dark lines). This is the phase diagram of three-dimensional 
supergravity. The tessellation by Tg (light lines) is shown to guide the eye. 
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Within this modular group Tg, there is a subgroup isomorphic to Z generated by 
the translation r — > r + 2. This Z is the intersection of Tg with the mapping class 
group of the handlebody Mo,i, so modular transformations in Z do not generate new 
handlebody contributions to the supergravity partition function. There is a unique saddle 
point contribution to F(r) for each element of the coset r#/Z. The subtessellation of the 
upper half plane Ji associated to this coset Tg/Z is shown in Fig. 4b). As one moves 
between tiles in this diagram, we expect a first order Hawking-Page phase transition as 
different saddles become dominant. So the subtessellation depicted in Fig. 4b) should be 
interpreted as the phase diagram of three dimensional supergravity. 

In |13|], holomorphically factorized partition functions for three-dimensional gravity 
were considered. Numerically, it appears that! the zeroes of these partition functions lie 



on the dark curves in Figure 4b) , and become dense on these curves in the large k limit |4(| . 
This indicates that the Hawking-Page transition in supergravity occurs by the mechanism 
of Lee- Yang and Fischer, as in the bosonic case described in Sec. 6. It may be possible to 
prove this by extending the analytic arguments of []37j and Sec. 6. 
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